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Abstract High resolution X-ray spectroscopy of the hot gas in galaxy clusters has 
shown that the gas is not cooling to low temperatures at the predicted rates of hundreds to 
thousands of solar masses per year. X-ray images have revealed giant cavities and shock 
fronts in the hot gas that provide a direct and relatively reliable means of measuring the 
energy injected into hot atmospheres by active galactic nuclei (AGN). Average radio jet 
powers are near those required to offset radiative losses and to suppress cooling in isolated 
giant elliptical galaxies, and in larger systems up to the richest galaxy clusters. This 
coincidence suggests that heating and cooling are coupled by feedback, which suppresses 
star formation and the growth of luminous galaxies. How jet energy is converted to 
heat and the degree to which other heating mechanisms are contributing, eg. thermal 
conduction, are not well understood. Outburst energies require substantial late growth of 
supermassive black holes. Unless all of the ~ 10^^ erg required to suppress star formation 
is deposited in the cooling regions of clusters, AGN outbursts must alter large-scale 
properties of the intracluster medium. 
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1 INTRODUCTION TO X-RAY CLUSTERS OF GALAXIES 



Clusters of galaxies are the largest gravitationally collapsed concentrations of 
matter in the Universe. With diameters of several megaparsecs and masses up 
to 10^5 M0, they are recognizable in photographs as distinct concentrations of 
galaxies centered on one or more brightest cluster members. The space between 
the galaxies is filled with a hot, dilute, plasma that emits X-rays and all is held in 
place by the gravity of a dark matter halo. The Perseus and Coma clusters were 
among the first clusters to be identified as X-ray sources by the Uhuru satellit e 
in the early 1970s ((Giacconi et al. 1971 : Gurskv et al. 1971 : Forman et al. 19721 ). 
By the mid-1970s, at least 40 clusters of galaxies w ere identified as extended 
and powerful X-ray sources (jCurskv fc Schwartd[l977l ). A thermal origin for the 
X-ray emission was confirmed in the Perseus, Virgo, and Coma clusters with the 
detectio n of the collisionally excited, 6-7 keV Fe-K emission feature by th e 
Ariel 5 (jMitcheh et al.lll976l ) and OSO-8 observatories (jSerlemitsos et al.lll977l ). 
Gas temperatures typically lie in the range of 10 million to 100 million K, cor- 
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responding to X-ray luminosities of Lx ~ 10^^ erg s~^ to 10^^ erg s~^. Roughly 
90% of the baryons in c lusters reside in the hot plasma, while the rest are locked 
up in stars in galaxies ( Lin et al. 20031 ). The hot gas forms a hydrostatic atmo- 



sphere, where the temperature and density distributions reflect the gravitating 
mass. The atmosphere serves as a bank of baryons that failed to end up in stars 
and galaxies, and as a repository for the heat exhaust and detritus from stellar 
evolution and the growth of supermassive black holes during galaxy and cluster 
formation. 

This review focuses on the latter, emphasizing new results from the Chandra 
and XMM-Newton X-ray Observatories showing that active galactic nuclei (AGN) 
lying at the hearts of galaxy clusters are pouring vast amounts of energy into the 
hot gas, some as prodigiously as quasars. The combination of high-resolution X- 
ray and radio imaging is yielding reliable measurements of this energy, which is 
apparently sufficient to suppress cooling flows and the substantial growth at late 
times of giant elliptical (gE) and cD galaxies. Deep Chandra images show that 
many clusters and gEs have multiple cavities, giving the hot atmospheres a Swiss- 
cheese-like topology that reveals the AGN outburst history independently of radio 
emission. These spectacular images are giving new insight into the particle and 
magnetic field content of radio sources, and are guiding the development of new 
radio jet and galaxy formation models. We discuss recent developments that link 
AGN outbursts to heating of the intracluster gas, and we tie these developments 
to several outstanding problems including the quenching of cooling flows, the 
exponential decline in the numbers of bright galaxies, the relationship between 
bulge mass and black hole mass in galaxies, and the possible contribution of AGN 
to excess entropy (preheating) in the hot atmospheres of groups and clusters. 
We briefly discuss some of the cosmological issues related to AGN heating in 
clu sters. General revi ews o f clusters f rom a n X-ray perspective were given recently 
bv iMushotzk^ (l2004l\ and jArnaudl Hooi), a review of clusters as cosmological 



probes was given bv IVoitj (Eo05l'l. and cold fronts and shocks associated with 



cluster mergers are reviewed by Markevitch &: VikhlininI ( 20071 ). We begin with a 



brief overview of the basic properties of clusters emphasizing the hot intracluster 
medium (ICM) and the scaling relations that describe it. 

1.1 X-ray Emission from the Intracluster Plasma 

The intracluster plasma (which we refer to interchangeably as hot gas) is com- 
posed primarily of ionized hydrogen and helium, mixed with traces of heavier 
elements, at roughly 1/3 of solar abundances. The presence of the gas can be 
understood in the context of hierarchical structure formation models. The warm 
baryons were swept inward with collapsing dark matter and subsequently heated 
to the virial temperature of the halo by accretion shocks and adiabatic compres- 
sion. Mean gas temperatures of several tens to one hundred million Kelvin reflect 
the virial temperatures of halos, so that T oc o"^, where a is the line-of-sight ve- 
locity dispersion of the cluster galaxies. Observed particle densities range from 
10~^ cm~'^ in the halos of clusters up to 10~^ cm~^ and higher in the centers of 
some clusters. 

The gas can generally be treated as an optically thin coronal plasma in ion- 
ization equilibrium. The electrons and ions interact through Coulomb collisions 
and radiate m ainly by thermal bremsstrahlung emission in the X-ray band (e.g., 
SarazinI 19881 ). At temperatures below 3 x 10^ K, X-ray emission from the gas is 
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increasingly dominat ed by th e recombination lines of iron, oxygen, silicon, and 
other elements (e.g., ISarazi 3 [1988). The emission from heavy elements, partic- 



ularly the iron K lines at 6 - 7 keV and the iron L lines below 1 keV in cooler 
plasmas, significantly alters the shape of the spectrum. This and the exponential 
decline in emission at high energies permit the temperature and metallicity of 
the hot gas to be measured accurately with modern X-ray telescopes. 

The X-ray telescopes used to study clusters have commonly employed propor- 
tional counters or charge coupled devices (CCDs) as detectors. They are sensitive 
to photons with energies spanning the range 0.1-10 keV, well-matched to ther- 
mal radiation from the intracluster gas {kT = 1 keV for T = 1.16 x 10^ K). 
Because the emission processes are collisional, the power radiated per unit vol- 
ume is proportional to the square of the density. X-ray surface brightness can 
therefore be used to determine gas density. With most X-ray instruments, for 
gas temperatures kT ~ 2 keV or greater, the count rate from an optically thin 
thermal plasma depends almost exclusively on its emission measure, J ngnn dV, 
where Ue is the electron number density, nn is the equivalent hydrogen number 
density, and the integral is taken throughout the emitting volume. As a result, 
gas densities can be determined quite accurately, even when gas temperatures 
are poorly known. 

The surface brightness of the hot gas declines with increasing radius, approxi- 
mately as Ix oc r^^ at large distances. Despite the rapid decline in surface bright- 
ness, the relatively low X-ray background permits X-ray emission to be traced 
to very large radii, making it an excellent probe of gas temperature, metallicity, 
and mass throughout much of the volume of a cluster. 

Surface brightness profiles have traditionally been characterized using the isother- 
mal "/^-profile" 

Ix « [1 + (r/rc)2]-3/5+i/2^ (1) 
where Tc is the core radius of the gas distributi on, and the parameter /3 ~ 2/3 for 



relaxed, bright clusters. As intr o duced (e.g.. ICavaliere &: Fusco-Femiand Il976l : 
Branduardi-Ravmont et al. 1981 : Forman JonesI 19821 ) . d is the ratio of the 



energy per unit mass in galaxies to that in the gas. However, this is an approx- 
imation that further relies on both the gas and dark matter being isothermal. 
This model provides a reasonably good fit to the surface brightness profiles of 
clusters at intermediate radii. In the central regions of some clusters, where the 
gas temperature declines and the density rises rapidly, the fit is poorer. At large 
radii, beyond roug hly O.Sr^nn, the obser ved surface brightness profiles steepen 



below the /3-profile ( Vikhlinin et al. 20061 ). Here, r2oo is the radius within which 



the mean mass density of the cluster exceeds 200 times the critical density of the 
Universe. For isothermal gas (or gas with kT > 2 keV), the /9-model corresponds 
to the electron density profile 

ne(r)=no[l + (r/rc)2]-=^^/^ (2) 

where no is the central electron density. 

More general ly, gas density and temperature profiles can be determined by 
"deprojection" ( Fabian et al. 198ll ). Under the often inadequate assumption of 



spherical symmetry, the X-ray spectrum at any point in the cluster is determined 
in terms of simple integrals of temperature and density. Typically, the gas is 
represented as a number of shells of uniform density, temperature, and compo- 
sition and its properties are de termined by fitting X-ray spectra, extra cted from 
corresponding annular regions (|Ettorill2000l : IPointecouteau et al.l[20o3 ). 
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1.2 Gas Temperatures and Masses 

The gas temperature of a hot atmosphere is most sensitive to the gravitating mass 
profile and, to a less er degree, the history of heatin g by gravitational and nongrav- 
itational processes (|Voitil2005l : iBabul et alJ 120021 ') . The connection between halo 
mass and gas temperature is clearly evi dent in the correlation between galaxy 
velocity dispersion and gas temperature ( Edge &: Stewart 199ll ). The scaling of 
this relationship is a oc T^-^^ in cluste rs with temperatures between 0.5 keV 
and 10 keV (e.s.. iKochanek et alJlioO^ ). which is close to the expected scaling 
a (X T^'^. Simple theoretical models of gravitational collapse predict the present 
day scaling relationships betwe en virial mass. X-r ay gas temperature, and lumi- 



nosity; M oc r3/2 and L oc (jEvrard et al.lllQQfil ). The observed scaling makes 



measured cluster temperatures and luminosities valued proxies for the much less 
accessible masses of dark matter halos. Departures from these scaling relations 
reflect physics beyond pure gravitational dynamics, including heating agents such 
as supernova explosion s and AGN, and additional pr essure support from magnetic 
fields and cosmic rays (|Markevitchlll998l : IVoitjl2005l ). 

When the gas is spherically symmetric, in hydrostatic equilibrium, and only 
the thermal gas pressure is significant, the run of temperature and density with 
radius, in principle, permits the mass profile of clusters to be computed from the 
equation of hydrostatic equilibrium as 



M(r) 



kTr 
Guniii 



dlogUe dlogT 
d log r d log r 



(3) 



where G is the gravitational constant, mn is the hydrogen mass, and fi ~ 0.61 
is the mean molecular weight of the ionized plasma. In practice, the logarithmic 
derivative of gas density is often evaluated using a parametrized model, such 
as the /5-model, for the gas density, with the parameters determined by fitting 
the surface brightness profile. Similar methods can be applied to finding the 
logarithmic derivative of the temperature. 

Mass determinations rely critically on the assumptions that the X-ray emitting 
gas is hydrostatic and the sole s ource of pressure. Other possible sources of pres- 
sure, including magnetic fields ( Govoni &: Feretti 20041 ) . cosmic rays, and bulk 
motion of the gas, will cause masses to be underestimated. Velocity dispersion 
and weak lensing masses generally agree with X-ray masses to within a few tens 
of percent, except perhaps in systems undergoing major mergers. Simulations 
have shown that bulk motions cause the h ydrostatic niass ap proximation to be 
bi ased below the true va lue by 5% to 20% (jNagai et al.ll2007l ). 

Vikhlinin et~aD (I200(tI ) used Chandra observations extending to large radii to 



measure the mass versus temperature relationship for relaxed clusters with tem- 
peratures in the range of 0.7-9 keV. They found M oc T^-^-i-^, which agrees 



with self-similar models ( Vikhlinin et al.l 20061 ). There are disagreements between 



researchers about the normalization constant and temperature profile shapes, re- 



Markevitch et al 



suiting in discrepancies in rnass by a few tens of percent fe.g 

19981 : llrwin fc Bregmanll2000l : lAllen et al.ll200ll : IPe Grandi fc Molendill2002l : IVikhhnin et al 
20061 ). These discrepancies are often associated with the measurement of dT/dr 



at large radii, where systematic effects dominate. However, some of the scat- 
ter may be real and perhaps associated with accretion, mergers, or AGN-related 
activity. 

The observed X-ray luminosity versus temperature relation for 3-10 keV 
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clusters scales as Lx oc r^-S-^-S which is steeper than expect ed for self-similar 
scaling (tx Markevitch et al. 19981 : Arnaud &: Evrard 19991 ). The steepening 
i s less significant in hot, relaxed clusters with gas temperatures above 5 keV 
(lAllen et alJl2nnih . )ut becomes pronounced in cooler clusters and groups. This 
departure has been attributed to excess entropy or preheating of the gas prior 
to virialization. Likely energy sources are the same supernovae that enriched the 
gas with metals and AGN outflows from nascent supermassive black holes. 

In the centers of clusters where radiati ve cooling is important, the gas temper - 
ature declines by factors of two to three ( Allen et al.ll2001 : Vikhlinin et al.r2005l ). 
The temperature there reflects the reduction in GM{r)/r over the central galaxy. 
It is in these so-called cooling flow regions where AGN outbursts are seen to have 
their biggest impact, as discussed below in detail. 



1.3 Mass Partitioning 

A galaxy cluster's gravitational influence extends over hundreds of cubic mega- 
parsecs, which is enough in principle to capture a representative sample of the 
dark matter and baryons in the Universe. The form of the matter and its dis- 
tribution throughout halos depends on both the history of galaxy formation and 
the cosmological parameters. X-ray and optical/IR observations show that the 
bulk of baryons reside in the hot ICM and that the fraction of baryons in stars 
decreases with cluster mass (the cold gas content is negligible). Only about 14% 
of the baryons in a 10^-^ M© cluster are in the stars, and this fraction decreases 
to 9% in 10^^ M0 cluster s, possibly owi ng to a declining efficiency of galaxy 
formation with halo mass (iLin et al.ll2"oo3 ). 

The gas mass fraction within rsoo, the radius within which the mean mass 
density of a cluster exceeds 500 times the critical density of the Universe, de- 
pends wea kly on gas temp erature, rising to ~ 12% in clusters hotter than ~ 4 
keV (.Vikhlinin et al.l l200fil ). Gas mass fraction also varies with radius, and is 
affected by any process that modifies the amount of star formation or the en- 
ergy content of the re maining gas, including radiative cooling and (pre-)heating 



(|Kravtsov et al.ll200Rl ). The baryon fraction in clusters is the sum of the gas 



fraction and the stellar baryon fraction. Based on near-infrared luminosities of 
galaxies in clusters with masses of ~ 3 x 10^^ M© (within rsoo), the mass fraction 
in stars is approximately 1.64%. Therefore, the cluster baryon fraction within 
this radius is ~ 14 %. This figure is clos e to the WMAP microwave background 3- 
year measurement ( Spergel et al. 200?! ). which gives a Universal baryon fraction 
of 16% to 19%. 

If the true baryon fraction is constant at a fixed radius in clusters and is 
constant over cosmic time (redshift), the dependence of the measured baryon 
fraction on luminosity distance can be a useful probe of the history of expansion 
of the Universe. Applying these assumptions to a study of 26 X-ray clusters 
extending to z = 0.9 with the Chandra observatory, lAllen et al.l ( 20041 ) found a 
high value for Ob that is consistent with type la supernova and WMAP3 values. 
It must be emphasized that the assumptions that this technique relies on have 
not been shown to be valid. The baryon fraction is affected by the history of 
radiative cooling, fee dback, and star forma tion, particularly in the cD, which are 
not well understood ( Kravtsov et al. 20051 ). 
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1.4 Magnetic Fields 

Faraday rotation measurements of background radio galaxies and radio galax- 
ies within c lusters have revealed that the intracluster gas is threaded with mag- 
netic fields (ICarilh k Tavloill2002l : [Kronberd[20o3 : lGovoni Feretti II2004I : IValled 
20041 ). The Coma cluster's ma gnetic field was one of the first to be detected at a 
level of a few microgauss fe.g.. lKronb"^l200.^ ). However, in the cores of cooling 
flow clusters (Section 2) field str engths of tens of rn icrogauss have been inferred 
using Faraday rotation measures (jClarke et alJl200lh . The ratio of magnetic pres- 



sure (-B^/Svr) to gas pressure {2nekT) is typically a few percent, so that magnetic 
fields are dynamically unimportant in clusters, except occasionally in the inner 
several kiloparsecs or so. Magnetic fields may have been deposited by radio galax- 
ies and or quasars, or they may be primordial fields that have been a mplified over 
time by gas turbulence and/or compression (ICarilh fc Tavloill2002l ). 



1.5 Transport Coefficients 

Magnetic fields can modify the transport coefficients significantly because the 
ratio of the Coulomb mean free path to the Larmor radius of a thermal proton is 
~ 10^5(A;T)3/2n-\ where the magnetic field, B, is in fjG, the temperature, kT, 
is in keV, and Uq is the electron density in cm~^. The same ratio for electrons 
is a factor of y^mp/m-c larger, so that electrons and ions are both tied rigidly 
to magnetic field lines for typical values of the intracluster magnetic field (e.g., 
GovonillioO^ ). One m ght then expect thermal conduction and vi scou s forces to 
be controlled entirely by the structure of the magnetic field (e.g., liable IQSi). 
However, under similar conditions in the solar wind, thermal conduction appears 
to be suppressed by only a modest factor an d it has been argued that this should 
be true generally for turbule nt plasmas (e.g., Sarazinlll988l : Rosner &: Tuckerlll989l : 
Naravan MedvedevlboOll ). 

Many empirical arguments favor greater suppression of thermal conduction in 
clusters. A number of researchers argue that thermal conduction must be sup- 
pressed from th e Braginskii value for a nonmagnetized plasma ( Spitzer 19621 : 
Braginskii 19651 ) b y factors of up to 1000 or more in order to h ave sharply de- 
fined "c old fronts" (jEttori fc Fabianll200d : iKemoner et al.ll2002l ). IVikhlinin et~al 



(|2001bl ) argue that cold fronts are a special case, because shear can amplify mag- 
netic field parallel to the front, effectively suppressing heat flow across the front. 
For example, they conclude that a ~ 10 //G field is required in the cold front 
in Abell 3667 to expla i n sup pression of Kelvin-Helmholtz instability, although 
Churazov &: Inogamovl ( 20041 ) point out that curvature suppresses growth of the 
instability, significantly weakening this c onclus ion. 

As another example, Ivikhlinin et all » found ~ 3 kpc remnants of the 
interstellar medium in NGC 4874 and NGC 4889 at the center of the Coma clus- 
ter. In order for these to survive thermal evaporation, they found that thermal 
conduction must be suppressed there by a factor of 30 - 100 from the Braginskii 
value. Again, these are special cases which rn ay not be representative of the gen- 
eral ICM. However, iMarkevitch et al.1 (|2003l ) argue that temperature variations 
in Abell 754 require thermal conductivity to be suppressed by at least an order 
of magnitude. Their result is determined for 100 kpc scale regions throughout 
the cluster that are not associated with any special structures. 

There are fewer limits on the viscosity. Based on the morphology of Ha fil- 
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aments near the northwestern ghost cavity in the Perseus Cluster, iFabian et al 
(|2003bl l argue that the flow there is laminar, hence that the viscosity is sup- 
pressed by a factor of no more than 15. T his interpretation is supported by 
numerical simulations ( Reynolds et al. 20051 ). With electrons and protons both 
tied so tightly to the magnetic field lines, it is reasonable to expect similar levels 
of suppression for thermal conductivity and viscosity. Taken together, the obser- 
vational limits suggest suppression of thermal conduction by somewhat more than 
expected from theory and observations of the solar wind and, perhaps, a similar 
level of suppression of viscosity. The constraints are showing signs of converging, 
but there is no consensus. Unsatisfactory though it is, transport coefficients in 
the hot ICM are not well determined. 



1.6 Metal Abundances 

The intracluster plasma is enriched w ith heavy elernents to an average metallic- 
ity of roughly 1/3 of the solar value ( Arnaud et al. 19921 '). which is roughly the 
universal average mean metal abundance ( Renzini 20041 ). Elements heavier than 
helium are created by stellar evolution, particularly by core collapse supernova 
explosions of massive stars (SNe II), and by thermonuclear detonations of accret- 
ing white dwarf stars (SNe la), which also play a significant role in their dispersal. 
Core collapse supernovae presumably seeded the Universe with metals during the 
early stages of galaxy formation, while SNe la, which are associated with the late 
stages of stellar evolution of less massive stars, dominate metal production over 
longer time spans. Clusters of galaxies are close to being "closed boxes" and 
thus retain the memory of metal enrichment through star formation and stellar 
evolution. SNe la are major producers of iron, whereas SNe II produce high 
yields of the alpha elements (Si, S, Ne, Mg). X-ray measurements of the relative 
abundance of the metallic species are thus able to constrain the history of star 
formation. 

The iron mass in the hot gas in clusters correlates with the luminosity of 
the elliptical and lenticular galaxy population but not with the luminosity of 
spirals. Moreover, alpha element abundan ces relative to iron in hot clusters are 
inconsistent with those in the Milky Way teaumgartner et al.l 2005 ). Therefore, 
most of the heavy el ements were created by stars bound to the early-type galaxies 
(lArnaud et al.lll992l ). The iron mas s in the gas exceeds the iron mass in stars in 
galaxies by at least a factor of two ( Arnaud et al. 19921 : Renzinil 20041 ). implying 
that cluster galaxies have ejected most of the metals they produced over cosmic 
time. If the iron in the ICM was shed by the stars, then the ratio of the iron 
mass to total stellar mass in clusters is a factor of four larger than expected for a 
population of stars like those in the Milky Way (|Lowensteinll200fil ). implying that 
star formation as it proceeds in the Milky Way would have difficulty producing 
the observed levels of iron in clusters. It is unlikely that SNe la could have 
supplied most of the iron unless the supernova rate was much higher in the past. 
Furthermore, SNe II underproduce iron by a factor of five if galaxies formed 
with a Salpeter initial mass function (IMF). It appears that most of the iron was 
produced by SNe II in a rapid phase of early star formation that proceeded with 
an IMF heavily weighted to ward massive s tars relative to th e IMF operating in 
present day spiral galaxies (Renzini 2004 : Lowenstein 20061 ). No evidence has 
been found for strong changes in abundance with lo okback time out to z ~ 0.5, 
which is consistent with early enrichment scenarios (jMushotzkv &: Loewenstein 
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19971 ). 

Abundance gradients corresponding to average metallicity increases of factors 
of two or more are routinely found in the central ^ 100 kpc regions of cD clusters 
(lEzawa e t allll99 7': 'ikebe et alJ|l997|; iFinoguenov et alJ[200ol : be Grandi Molendi 
I2OOI : Irwin Sz Bregman 200L ). High-resolution Chandra studies have shown that 
the gradients are strongest nea r the cD and that they someti mes track the stel - 
lar isophotes ( Wise et al. 20041 ) . For example, in Hydra A, ( David et alJ 200ll ) 
found that the iron abundance rises from 30% of the solar value at 100 kpc to 
about 60% in the nucleus of the cD, and that the silicon abundance rises from 
roughly half the solar value to the solar value in the nucleus. They attribute 
the gradient to constant SNe la production from the cer itral galaxy over the pas t 
several gigayears. A similar situation is found in M87 (jFinoguenov et al.ll2002l '). 
where most of the iron seems to have originated from SNe la, while only 10% 
was injected from SNe II. Intriguing ly, iTamura et aP (|2004l ) found that SNe II 
enrichment from massive star formation may be important in the cores of some 
clusters. Star formation is often observed in cD galaxies centered in cooling flows, 
which could contribute significantly to the gradients. 



2 CLASSICAL COOLING FLOWS 



A cooling flow cluster is characterized by bright X-ray emission from cool, dense 
gas in the central region of the cluster (see Fabian 19941 for a more comprehensive 
review of the basic operating principles of a cooling flow). Within the cooling 
radius, where the cooling time of the gas is less than the time since the last major 
heating event, the surface brightness of the gas near a central cD galaxy often 
rises dramatically above a /3-model, by factors of up to 100, corresponding to a 
rise in gas density by factors of 10 or more. The X-ray luminosity within the cool- 
ing region reaches values of 10^^ erg s~^ in the extreme. In many cases it is more 
than 10% of the cluster's total luminosity. If this luminosity is uncompensated 
by heating, the gas will radiate away its thermal an d gravitat i onal energy on a 
timescale of tmn] = p/\H — l)ne,rmA(T)] < 10^ yea r (ISiSll976l : rCowie BmnevI 
I977I : iFabian fc Nulsenll 19771 : iMathews fc BregmanllT97i r where p is the gas pres- 
sure, A(T) is the cooling function, and 7 is the ratio of specific heats of the gas. 
As the gas radiates, its entropy decreases and it is compressed by the surround- 
ing gas, causing it to fiow inward. The cooling time decreases as the gas density 
increases and, eventually, the gas temperature drops rapidly to < 10^ K, so that 
cooled gas condenses onto the central galaxy. The condensing gas is replenished 
by hot gas lying above, leading to a steady, l ong-lived, pr essure-driven inward 
flow of gas at a rate of up to 1000 Mq year" -1 (|Fabianlll99^ 1. 

Observed cooling times are significantly longer than free-fall times, so that the 
gas remains very nearly hydrostatic as it cools. The fiow is then governed by 
cooling, making the flow time t^o^ = r/v, where v is the radial speed of inflow, 
approximately equal to the cooling time. Counterintuitively, the heat lost to ra- 
diation does not necessarily make the gas temperature decrease. The inexorable 
entropy decrease is offset by adiabatic compression as gas flows inward, so that, 
typically, the temperature of the cooling gas follows the underlying gravitational 
potential, i.e., kT/(fim}i) is a multiple of order unit y of the "local virial tem- 



perature," GM{r)/r (e.g., solutions of lNulsenlll986l with k = 0). Cooling gas 



would therefore remain approximately isothermal in an underlying isothermal 
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potential. Because the virial temperature of cluster central galaxies (revealed by 
their stellar velocity dispersions) is typically lower than that of the surrounding 
cluster, the gas temperature declines inward in a classical cooling flow. Rather 
than being a direct manifestation of cooling, this temperature drop reveals the 
flattening of the underlying gravitational potential. (If gas cools enough to flow 
inward, this feature will persist in more up-to-date cooling flow models.) The 
Mach number of the flow increases inward, until it approaches unity. Up to that 
point, linear growth o f the thermal instability is suppressed by buoyant motions 
( Balbus &: Soker 19891 ). but beyond it cooling is too fast and the gas is expected 



to form a rain of thermally unstable clouds that cool rapidly to low temperature. 

The classical cooling flow is approximately steady within the region where the 
cooling time is shorter than its age. The power radiated from the steady flow 
equals the sum of the enthalpy carried into it and the gravitational energy dissi- 
pated within it. To a first approximation, the gravitational energy can be ignored, 
so that the luminosity Lx ^ M{5kT)/{2iJ,mn) ^ 1.3 x 10'^^r5M22 erg s'^, where 
the temperature of the gas entering the fiow is given by T5 = {kT/5 keV) and 
the cooling rate is given by M2 = M/(100 Mq year"-*^). This expression is exact 
for steady, isobaric cooling. 

This single-phase cooling model predicts central spikes in surface brightness 
that are clearly stronger than observed, a discrepancy that prompted the intro- 
duction of inhomogenous cooling flow models. Such mod els postulate a broad 
spectrum of gas temperature and density at each radius (Nulsen 19861 ). Both 



Rayleigh- Taylor and shear instabilities can disrupt an overdense cloud if the dis- 
tance it must move to reach its convective equilibrium position exceeds its size 
(alternatively, if its fractional overdensity is greater than its size divided by the 
distance to the cluster center). Such nonlinear overdensities therefore tend to 
be shredded finely, slowing their motion relative to the bulk fiow and enabling 
them to be pinned to the flow by small stresses (e.g., magnetic stresses). Tying 
overdensities to the bulk flow then permits the growth of thermal instabilities 
throughout the core of the cluster, leading to widespread deposition of cooled 
parcels of gas and making the mass flow rate, M, a strong function of the ra- 
dius. This solution avoids the need for strong X-ray and UV surface brightness 
spikes centered on the cD's nucleus by distributing the hundreds of solar masses 
of gas and star formation deposited each year throughout the cooling region of 
the cluster. 

This cooling flow model has been under siege for years, primarily because of 
its failure to predict the observed amount and spatial distribution of star forma- 
tion, line emission, and other expected products of cooling, which are generally 
observed only in the inner few tens of kiloparsecs. The spatial distribution of 
cool gas and star formation is more consistent with the single-phase model, but 
at levels that fall orders of magnitude below the predictions. This failure implies 
that the gas is not condensing at the predicted rates, and that radiation losses are 
either being replenished, or the gas is condensing into an unseen state. Sensitive 
searches for the repository in optical, infrared, and radio bands have severely 
restricted the latter possibility, if not ruled it out entirely. 

2.1 The Modern View of Cooling Flows 

The strongest spectroscopic signatures of co oling gas are X-ray em ission lines 
below 1 keV of various charge states of Fe L (jBohringer et al.ll2002l ). Early ob- 
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servations of M87 and a few other bright clusters m ade with the Eiristein Obser- 
vatory's Focal Plane Crystal Spectrometer [FPCS] (jCanizares et al.lll98l Il988l ) 
apparently detected these Fe L lines at the 4cr to 6a level. Remarkably, the line 
strengths from gas at temperatures of a few million to a few tens of million de- 
grees agreed with predictions of the classical cooling flow model. Unfortunately, 
the FPCS, with an effective area of only one square centimeter, was unable to 
achieve detections of higher statistical significance. Nevertheless, for two decades 
these concordant results were bolstered by signatures of c ooling at other wave- 
lengths: fil aments of Ha ernissioii from warm, ionized gas (Heckman 198ll ): star 



formation (jJohnstone et al.l Il987l: iMcNarnara &: O'Connelll Il989l ): and pools of 
cold atomic and molecular gas "(|Edgell200ll ')" provided support for a classical cool- 
ing flow model, albeit at lower rates. 

The situation changed dramatically when far more sensitive XMM-Newton 
Reflection Grating Spectrometer (RGS) observations failed to conflrm the picket 
fence of lines, including the Fe L features, at ^ 1 keV from cooling gas in spectra 
of th e clusters Abell 1835 (jPeterson et al.1 12OOII ) and Abeh 1795 (jXamura et al. 
200 ll ). The prominen t and useful Fe XVII line at 12 A is either weak or absent in 
cooling flow clusters ( Peterson et al. 20031 ) . although there is a hint of Fe XVII 



emission in the Abell 2597 cluster at a level that is consistent with a cooling rate 



of ~ 100 Mq year" -1 (|Morris &: FabianI l2005l l . Ultraviolet line emission also sug- 



gests cooling rates about an order of mag nitude smaller than previous estimates 



(jOegerle et al.ll200ll : iBregman et al.ll200' 



The relationship between cooling rate and line power that is central to this issue 
was flrst noted by Cowie (1981). In terms of the entropy, S, the energy equation 
of cooling gas is pTdS/dt = —ncniiA{T), where p is the gas density, t is the 
time and other terms are as above. Rearranging this to describe a fluid element 
in gas cooling steadily at the rate M gives the mass of the element as dM = 
M dt = —MpTdS/[neniiA] = pdV, where dV is its volume, so that the emission 
measure of the element is rienn dV = —MTdS/A (where dS < for cooling gas). 
In terms of the pressure and temperature, dS = [7 dT/ (7 — 1)T — dp/p]k/ (pmn). 
As the gas cools and flows inward, its pressure generally increases, so that —dS 
is minimized for isobaric cooling (dp = and dT < 0; late stages of cooling can 
be isochoric, reducing the integrand below by a factor of 7 at sufficiently low 
temperatures). If the contribution of a line (or group of lines) to the cooling 
function is Ai[^c(T), then the total power radiated in that line by gas cooling 
isobarically in a steady flow is 



line 



nenn Aline (T) dV = M 



7 



k 



7 - 1 /iTTlH 



A(r) 



dT, 



where T is the temperature that the gas cools from (the sense of integration has 
been reversed here). If the gas pressure increases as the gas cools, the power 
radiated in the line would be greater and this expression would overestimate the 
cooling rate. 

The failure to detect the low energy X-ray lines at the expected levels indicates 
that the canonical cooling rates were overestir nated by an order of mag nitude or 
more, although the picture is not so simple ( Peterson Sz Fabian! 120061 ). In the 
classical cooling flow model, thermally unstable gas clouds are expected to be 
approximately isobaric, so that the constant pressure cooling model should be 
reasonably accurate, especially for lines emitted mainly by gas that is cooler 
than average. As outlined above, the constant pressure cooling model makes spe- 
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cific predictions of the strengths of lines from low temperature cooling gas (cf . 
Bohringer et al. 20021 : Peterson &: Fabian 2006 ) . These predi ctions are inconsis- 



tent with the line strengths observed by Peterson et al. ( 20031 ) . which require the 



emission measure to decrease more rapidly with temperature than expected for 
the constant pressure cooling flow model. This behavior cannot be explained by 
merely reducing the cooling rate. 

The absence of multiphase gas is confirmed in moderate spectral resolution 
CCD data from Chandra, XMM-Newton, and ASCA. In general, the deprojected 
gas profiles in cooling flows can be adequate l y modeled by a single temperature 
plasm a at each radius ([McNamara et al. I l2000l :l lOavid et al.ll200ll : lMolendi &: Pizzolatol 
200ll ). except perhaps near the nucleus of the cD, where star formation is fre- 
quently observed. 

It may be possible for nonradiative cooling to proceed at or near the canonical 
rates without revealing X-ray line emission. Processes such as gas-phase mix- 
ing, differential photoelectric absorp tion, and rapidly c ooling, unresolved, high 
metallicity inclusions in the hot gas ( Fabian et al. 200ll ) are possible. If so, the 
problem of identifying the permanent repository for the cooling material would 
remain. 



RADIO LOBE RELATED X-RAY STRUCTURE AND CAV- 
ITIES IN CLUSTER CORES 



Disturbances in the hot gas near NGC 1275 were first noted in an early Ein 



stein Observatory i mage of the Perseus cluster (jBranduardi-Raymont et al.lll981 



Fabian et ai]ll98lh . A decade later, Rosat's 5 arcsec High Resolution Imager 



(HRI) associated the disturbances with tw o cavities filled with radio emission 



emanating from the nucleus of NGC 1275 (jBohringer et al.lll993l ). Similar dis 



turbances were later noted in HRI images of other bright cD clusters (e 
Carilh et al1ll994l : iHuang fc Sarazinlll998l : lOwen fc Eiie^ll998l : iRizza et al.lto 
but limitations in Rosat's spatial and spectral resolution stalled any further 
progress on the nature and import of these disturbances until the launches of 
Chandra and XMM-Newton in 1999. 

We now know that nearly three dozen cD clusters and a similar number of gE 



galaxies and groups harbor cavities or bubbles in their X-ray halos ( Fabian et aL. 

200d : iMcNamara et al.ll2000l . I2OOII : ISchindler et"al]l200ll : iHeinz et al]l2002l : lMazzotta et al 



20021 ) . Cavity systems are difficult to detect, so this is surely a lower limit to their 
numbers. Like the radio lobes that created them, cavities are usually found in 
pairs of approximately elliptical X-ray surface brightness depressions, 20% to 40% 
below the level of the surrounding gas. This is the expected decrement toward a 
spheroidal empty cavity embedded in a /3-m o del at i nosphere ( Fabian e t al.l 2000l : 
McNamara et al.l I2OO0I : iBlanton et aD I2OOII . l2003l : iNulsen et al.. ,2002 1. Cavity 
systems in cluster s vary enormously i n size , from diameters smaller than 1 kpc 
like those in M87 (" Forman et al.l2005 . 20071 ) to diameters approach i ng 200 kpc in 



the MS0735. 6+7421 and Hydra A clusters ([McNamara et al.ll2005l : iNulsen et al. 



2005b : Wise et al. II2007I ). A correlation exists between radio luininosi ty and cav- 
ity power (Section 5.1; e.g., Birzan et al. 2004 : Dunn &: FabianI 20061 ). but with 
a large scatter that is poorly understood. 

One of Chandra's early surprises was the discovery of cavities devoid of bright 
1.4 GHz radio emission. They were dubbed ghost cavities, and were interpreted 
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as aging radio relics that had broken free from the .jets and had risen 20 — 30 kp c 



into the atmosphere of the cluster ( McNamara et al. 2001 : Fabian et al. 20021 ) 



We now know they are filled with low-f requency radio ernission and may b e 
connected by tunnels back to the nucleus ( Clarke et al. 20051 : Wise et al. I 2007 ). 
Some clusters have multiple pairs of cavities, apparently produced by multiple 
radio outbursts or quasicontinuous outfiows. 

The work required to infiate the ca vities against the surrouri ding pressure 
is roughly pV ~ 10^^ erg in gEs (e.g., 'Finogue nov &: Joned I2001I ) and upward 
of pV = 10^^ erg in rich clusters (e.g., Raffer tv et al.l 20061 ). ^The total en- 
ergy needed to create a cavity is the sum of its internal (thermal) energy, E, 
and the work required to infiate it, i.e., its enthalpy, H = E + pV. This is 
several times pV. The displaced gas m ass is several 10^° U p, m an average 



200 ih but can exceed 



cluster system such as Abell 2052 (e.g., iBlanton et al. 
10^^ Mq in powerful outbursts such as those in MS0735. 6-1-7421 and Hydra A. 



The cavities in these systems occupy between 5% and 10% o f the volume withi n 



300 kpc giving the hot ICM a Swiss cheese-like topology (jWise et al. I l2007l ). 
The brigh t rims or shells surrounding many cavitie s are cooler than the am- 
bient gas (iFabian et al.]l200d : iMcNamara et all I2OO0I : iBlanton et al.l[200ll . l2003l : 
iNulsen et al.ll2002l) and thus are not activ e shocks as anticipated in early models 



(e.g., iHeinz. Revnolds. &: BegelmanI Il998l ). Evidently, the cavities are close to 
being in pressure balance with the surrounding gas. A nearly empty cavity will 
rise into the cluster atmosphere like a buoyant weather balloon, traveling at a 
speed approaching the local free-fall velocity. The cool rims are probably com- 
posed of displaced gas dragged outward from the center b^ the buoyant cavities 



(|Blanton et al.ll200ll : IChurazov et al.ll200ll : iR.evnolds et al. 2001'). 



Ca yity systems are often surro unded by belts (ISmith et al.„2002l). arm s (lYoung et al 



2OO2I : iForman et al.ll200,5l . l2007l ). filaments, sheets (jFabian et al.ll2006l ). and frag 



ile tendrils of gas maintained agains t thermal evaporation, perhaps, by magneti c 



fields threaded along their lengths ( Nipoti &: Binney 2004 : For man et al. 20071 ) 



This structure is usually composed of cooler gas and is associated with H emis- 
sion that may be tr acing circulation driven by rising radio lobes and cavities 



(jFabian et al.l l2003bl ). ' 'Swirls" of cool X -ray gas were foun d in the central re 



gions of the Perseus (jFabian et al.1 12OO6I ) and Abell 2029 (jClarke et al.l l2004l ) 



clusters, which may be related to merger activity or circulation fio ws generated 
by radio sources (e.g.. lMathews fc Brighelrtill2003l : iHeinz et al.ll2006l ). At present, 
this wealth of structure is not well understood. 

Owing to their proximity and high surface brightnesses, the M87 and Perseus 
clusters are spectacular exemplars of these structures (Figured]). Keep in mind, 
however, that neither is outstanding in its power output. The AGN (cavity) 
power in Perseus is about 25 times larger than that of M87, and Perseus itself 
is feeble compared to Cygnus A and MS0735. 6-1-7421, w hich are roughly 21 5 



and 270 times, respectively, more energetic th an Perseus (Raffertv et al. 20061 ). 



Some 18 of the 33 cavity systems studied bv iRaffertv et al.l ( 20061 ) exceed the 
AGN power output of Perseus, suggesting the Perseus and Virgo clusters serve 
as useful benchmarks for average to low power outbursts. 



3.1 Shock Fronts 



With the notable exceptions of the Mac h 8 shock in Centaur us A (jKraft et al. 
20031 ). the Mach 4 shock in NGC 3801 (jCroston et al.ll2007l ). and perhas sur- 



14 



McNamara &: Nulsen 



rounding the overpressured quasar 3C186 (|Siemiginowska et alJ [2OO5I I there is 
Httle evidence for strong shocks surrounding the radio-lobe cavities studied with 
Chandra. Howev er, evidence of weaker, more remote shocks is accumulating. 
Jones et al.l (120021) repo rted gas features near NGC 4636 consistent with a Mach 



1.7 shock. Smith et al. I 11002) saw evidence of shocks around the radio source 
Cygnus A. More recently, Chandra images of clusters have revealed roughly 
elliptical breaks in surface brightness enveloping the inner radio-lobes of M87 
(jForman et al.l I2OO5I ) and the en ormous cavity systems in MS0735+7421 ( Fig- 
ure[2]). Hydra A, and Hercules A (jMcNamara et al.|[2005l : iNulsen et al.ll2005a .b). 
These structures resemble the classical cocoon shocks thought to envelope pow- 
erful radio sources as they advance into the surrounding medium. Modeled as a 
spherical shock from a nuclear outburst, the surface brightness discontinuities can 
be reproduc ed by conventional sho c k waves with Mach numbers lying between 
1.2 and 1.7 (jMcNamara et al.ll2005l : iNulsen et"alll2005al lbl). Though the shocks 
are only mildly supersonic (weak) they can encompass huge volumes, 200 — 400 
kpc in radius and require energy deposition upward of 10^^ erg into the intraclus- 
ter gas. The shock energy corresponds to several pV per cavity, com parable to 
the free energy of the cavities themselves (see also I Wise et al. 1 120071 ) . Thus the 
cavity enthalpy provides only a lower limit to the total outburst power. The ages 
of these outbursts, estimated from shock models (approximately the shock radius 
divided by the shock speed), are reassuringly similar to buoyant rise times of 
their cavities, although significant differences found, for example, in Hydra A are 
related to the outburst history and the partitioning of energy between enthalpy 
and shock heating (Section 3.5). 

Detecting the expected gas temperature rises behind these shock fronts has 
proved difficult. Radiative losses are negligible on timescales of interest, so that 
shocks in the ICM should conserve energy. Furthermore, the magnetic field in 
the bulk of the ICM is dynamically insignificant and does not affect sho cks sig- 
nificantly for pure magnetohydro dynamic shocks (e.g., Nulsen et al.l20o3 ). Some 
shock energy can be absorbed by particle acceleration, but this is expected to be 
small for weak shocks. Under these conditions, the ICM should be well approxi- 
mated as an ideal gas with a constant ratio of specific heats, 7 = 5/3, leading to 
the well-defined relationship. 



T2 ^ (7 + l)P2/pi - (7 - 1) 
Ti [{■y + l)-{j-l)p2/pi]p2/pi 



(4) 



between the shock temperature jump, T2/T1, and density jump, P2I P\- For ex- 
ample, for the Mach 1.2 shocks in Hydra A and M87 the temperature jumps by 
20% behind the shock. However, the temperature declines rapidly behind the 
shock because of adiabatic expansion. When projected onto the sky, emission 
from the shocked gas is also diluted by emission from the surrounding unshocked 
gas. As a result, the projected emission-weighted temperature has a peak rise 
of only 5%. Given that thousands of photons are required to measure temper- 
atures accurately, few such temperature ri ses have been detec ted. Examples of 
detections include the 1 4 kpc ring in M87 (jForman et al.ll2007l ). MS0735.6-h7421 
( McNamara et all 20051 ) . Hercules A (PEJ Nulsen, in prepa ration), Centaurus A 
(jKraft et al.ll2003l l. and NGC 4552 (jMachachek et al.ll2006l ). 

For MS0735. 6-1-7421 and Hercules A, the mean jet power released into the ICM 
over ~ 108 years is more than 10^^ erg s ^, which is comparable to a powerful 
quasar. The bulk of this energy is deposited beyond the cooling region of those 
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clusters. Despite this quasar-like power, the cD hosts share few quasar character- 
istics (e.g., broad nuclear emission lines), although this issue needs further study. 
The energy deposited in the inner 1 Mpc of MS0735. 6+7421 corresponds to a 
few tenths of a keV per particle, which is a significant fraction of the ~ 1 keV 
required to supply the excess entropy (preheatin g) in clusters ( Wu et al. 2000l : 
Rovchowdhurv et al.l[2nn3 : IVoit fc Donahuell2nn,^ ) . Thus several outbursts of this 



magnitude during the life of a cluster, particularly in the early stages of its devel- 
opment, could preheat it (see Section 1.2). No clear distinguishing characteristics 
have yet been noted between the cD galaxies hosting powerful outbursts and cDs 
in other cooling flows. Thus occasional powerful outbursts, if they occur in all 
systems, could dominate the energy output from smaller, more frequent outbursts 
in tegrated over clu s ter ag es . 

Voit &: Donahuj ( 20051 ) have identified several clusters with strongly boosted 



central entropy profiles and relatively weak radio sources that otherwise show 
no evidence of recent AGN activity. They suggested that the central entropy 
boosts were imprinted by shocks generated by powerful AGN outbursts that 
occurred in the past, although the AGN are dormant at present. Furthermore, 
Voit & Donahue found that the entropy profiles are generally consistent with 
shock heating in the inner few tens of kpc, but the mode of heating switches 
to cavities and sound waves further out. An additional consequence of repeated 
powerful outburts is the substantial growth at late times of the black holes at the 
core of the AGN, which is discussed further in Section 7.2. 

Luminosity boosts, gas clumping, and gas outfiows caused by AGN outbursts 
could, in principle, affect measurements of cosmologically important quantities 
such as gas mass fractions, the luminosity function, and cluster masses. However, 
these effects are like l y to be subtle an d are just beginning to be explored (e.g., 



Kravtsov et"aDl2005l : [Gitti et al.ll2007l ) 



3.2 Ripples and Sound Waves 



The best known example of weak shocks or sound waves is that in the Perseus 
cluster. It is seen as a sp ectacular system o f ripples and other disturbances in an 
890 ksec Chandra image (|Fabian et al.ll2006l ). The half dozen ripples lying beyond 
the inner cavity system (Fi gure [J) are separated by roughly 11 kpc and are visible 



out to a radius of 50 kpc (jFabian et al 



l2003al ). The ripples appear to be pres- 



sure disturbances with amplitudes of 5% to 10%, or sound waves (weak s hocks ) 
propagating outward with a period of ~ 10 year (iFabian et al.l l2003al . l200fil ). 
Apparently the gas within 25 kpc surrounding the cavity system is overpressured 
by about 30%, implying that the entire inner halo is expanding. Cavity pV work 
alone would then underestimate the AGN power in Perseus. Temperature jumps 
across the ripples and the band of high pressure have not be en found, desp i te the 



Fabian et al 



(|2006l ) 



ample number (70 million) of detected photons available, 
have suggested that thermal conduction is suppressing the temperature jumps 
by creating isothermal shocks. This phenomenon is poorly understood. 



3.3 Detectability of Cavity Systems 

Studies of cavity populations should in principle yield information on the AGN 
duty cycle, the energy per AGN outburst, and outburst ages, once a reliable 
dynamical model has been established and the complicating issues of cavity de- 
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tectability, stability, and lifetime are understood. In a genuinely random sample 
of galaxy clusters, radio jets and their associated cavities are expected to emerge 
from the cD's nucleus at random orientations with respect to the plane of the 
sky. The decrement in surface brightness of a cavity relative to the surrounding 
cluster, i.e., its detectability, will be a strong function of its size and distance 
from the cluster center, as well as its aspect with respect to the plane of the sky. 
In the simplest approximation, a small bubble of radius r on the plane of the sky 
at a distance R from the cluster center produces a count deficit (cavity) scaling 
as r^(l + R? / a?)~'^l^ , where j3 and a are the /3-model parameters for the cluster. 
The cluster count from over the bubble scales as r^(l + R? / a?)"^^^^/"^ and the 
noise in this scales as its square root. Thus the signal-to-noise ratio scales as 



\l + R'/a') 



-3/3/2-1/4 



implying that cavities are easiest to detect when they 



are large and are located in the bright central regions of the cluster. This sim- 
ple argument does not take projection into account, but should apply roughly 
to bubbles within 45° of the plane of the sky as seen from the cluster center. 
Bubbles far from the plane of the sky are difficult to detect at any radius. 

The detectability of a cavit y system as a function of its age and nuclear dis- 
tance (time) was modeled by Enfilin &: Heinj ( 20021 ). who considered spherical 
bubbles of adiabatic, relativistic plasma embedded in an isothermal cluster with 
a core radius = 20 kpc. As a cavity rises, it expands adiabatically from its 
initial volume, Vq, to a modestly larger volume, Fi, at its observed location as 
= yo{Po/PiY^^ 1 where po and pi are the ambient pressures at the respective 
locations, and T is the ratio of specific heats within the cavity. For a spheri- 
cal cavity, the ratio of its radius, r, to the distance, R, from the cluster center 
therefore evolves as 



(5) 



where p{R) is the pressure of the ICM at distance R and subscripts "0" denote 
initial values. EnBlin & Heinz found that the contrast of a cavity launched from 
the center of the cluster along the plane of the sky is a slowly decreasing linear 
function of distance, until the cavity vanishes into the background at large dis- 
tances. The detectability of cavities rising at oblique angles with respect to the 
plan e of the sky dec l ines a ccording to a power law with distance. From figure 
3 of Enfilin &: Heinj ( 2002 ). we find that the detectability declines oc R^^ and 
oc R~^'^ for cavities launched at 45° and 90°, respectively, from the plane of the 
sky. These estimates do not include the efi^ects of cavity rims, which enhance 
cavity contrast, nor of cavity disruption (Section 3.6), that works in the opposite 
direction. Nevertheless, they imply that once cavities have ventured distances of 
several times their initial radii, their chances of detection decline rapidly. 
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3.4 Cavity Statistics 

The statistical properties of cluster cavities drawn from the Chandra archive, to- 
gether with their attendant radio sources , have been studied in some detail (e.g., 
Birzan et al]l2004l : lOunn fc FabianI l2004l . l2006l : lOunn et al.ll2005l : iRaflFertv et ah 
20061 ^. These studies are hobbled by unknown selection effects in the Chan- 
dra archive. Nevertheless, they suggest a high incidence of detectable cavity 
systems in clusters, groups, and galaxies, spanning a lar ge range of gas tern pera- 
ture. The detection frequencies found in the samples of Birzan et al. ( 20041 ) and 
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Raffertv et aD (|2006l l are 20% (16/80) and 25% (33/130), respectively. Using 



the flux -limited sarnple o f the brightest 55 clu sters observed by the Rosat obser- 
vatory ( Peres et alJll99^ ). lOunn et al. 1 (|2005l l found between 12 and 15 clusters 



with cavity systems, giving an overall detection rate of between 22 % and 27%. 



Thei r detection rate rose to 70% (12/17) in strong cooling flows (iDunn et al 
200.51 ). The detection rate is clo se to 25% (27/109) for a sample of nearby gEs 
with significant hot atmospheres ( Nulsen et al.ll2007l ). While there is considerable 



overlap between the three cluster samples, the gE sample is largely independent 

(Secti on 4.3). 

The Raffertv et al. ( 20061 ) sample is the most extensive cluster sample available. 



and thus provides a good representation of average cavity properties in clusters. 
The distribution of projected distances between the nucleus of the host cD and 
cavity centroid in Figure [3] shows that the detection rate peaks in the inner 30 
kpc or so and declines rapidly at larger distances. Only the rarest and most 
powerful outbursts produce detectable cavities beyond ~ 100 kpc. Within 100 
kpc the detection frequency declines formally as ~ R~^'^, but is consistent with 
~ R~^, the expected rate of decline for cavities launched on random trajectories 
(Section 3.3). 

The distribution of observed cavity sizes (Figure H]) shows a typical value of 
10 - 15 kpc in radius with no preferred size. The distribution of the ratio of 
projected nuclear distance to radius, shown in Figure O peaks at R/r ~ 2 and 
falls off rapidly beyond. Evidently, cavities travel roughly their own diameters 
before they disintegrate or become too difficult to detect. Other than the most 
powerful systems, cavities are found within the light of the central galaxy. The 
distribution of buoyancy ages (Figure [6]) shows a typical age of ~ 10'' years but 
some are greater than ~ 10^ year. 

The difficulty of detecting cavities coupled with the highly variable image qual- 
ity (depth) of the Chandra archive suggest that the existing inventory of cavities 
is incomplete. Most, if not all, cooling flow clusters probably harbor cavity sys- 
tems or have done so in the recent past. 

3.5 Cavity Kinematics and Ages 

During its initial stages, the tip of a radio jet advances supersonically into the sur- 
roun d ing medium (e.g., Scheu er'l974': Gu ll &: Northov er 1973: Bl andford &: Reed 



19741 : iHeinz. Reynolds, fc Begelman ,19981 : lEnBlin k Heinz ,2002 ). As the ram 



pressure of the jet declines with respect to the ambient pressure, it decelerates, 
eventually shuts down, an d the radio lobe s quickly reach pressure balance with 



the surrounding hot ICM (iBegelmanI l2004l ) . During the initial, supersonic stage. 



the cavity created by the tip of a jet can be long and narrow. Only ICM near 
the advancing jet tip is subjected to the strongest shocks, and relatively little of 
the X-ray emitting gas within the region encompassed by the jet is shocked. A 
narrow cavity produces a small deficit in the X-ray emission, making it difficult 
to detect. Enhanced X-ray emission from the relatively small volume of shocked 
gas may not stand out when superimposed on the general cluster emission. In 
order to form the roughly spherical cavities that are observed, the expansion of 
a cavity has to "catch up" with the jet tip, so that the cavity displaces a much 
greater proportion of the X-ray emitting gas within the extent of the jets. This 
may come about when the pressure within a radio lobe becomes comparable to 
the ram pressure of the jet, so that the lobe expands into the surrounding gas at 
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a speed similar to the t ip of the jet, or b ecause the jet wanders enough to carve 
out a large cavity (e.g., Heinz et al. 20061 ). At formation, such a cavity can have 
an age that is appreciably shorter than its sound crossing time and a lot less than 
its buoyant rise time. 

A cavity is buoyed outward with a force, -Fb = Vg{pa. — Ph), where V is its 
volume, g is the local gravitational acceleration, pa is the ambient gas density, 
and /9b is the density of the cavity. As the influence of the jet wanes, the cav- 
ity will start to rise outward under control of the bouyant force. Its terminal 
speed, determined by balancing the buoyant force against the drag force, is then 
vt ~ ^j2gV I (SC) ~ (4:V K/3)-y/2r/R, where S is the bubble's cross section, and 



C is the drag coefficient ( Churazov et al. 200 ll ). The second form is for a spher 



i cal cavity of radius, r, at a distance R from the cluster center, with C = 0.75 
( Churazov et al. 200 ll ). The Kepler speed, vk = y/gR, is comparable to the local 
sound speed, so that the terminal speed is almost invariably subsonic. In practice, 
the volume, V, is determined from X-ray measurements of the projected size of 
a cavity. Because the cavities generally lie within the body of the central galaxy, 
g can be estimated using the local stellar velocity dispersion, under the approx- 
imation that the galaxy is an isothermal sphere, as 5 ~ 2a'^ /R. Alternatively, 
the gravitating mass distribution can be determined on the assumption that the 
surrounding gas is hydrostatic, as outlined in Section 1.2, then used to calculate 
the gravitational acceleration. 

Three estimates are commonly used for cavity ages: the buoyant rise time, 
the refill time, and the sound crossing time. From above, the time taken for a 
bubble to rise at its buoyant terminal speed from the center of the cluster to 
its present location, i.e., its buoyant rise time, is approximately tbuoy — R/'Vt — 
RyJ SC / {2gV). This is a reasonable estimate for the age of a cavity at late 
times, long after it has detached from the jet that formed it. The "refill time" 
is the time required for gas to refill the displaced volume of the cavity as the 
bubble rises, i.e., the time taken for a cavity to rise buoyantly through its own 
diameter. If a cavity is formed rapidly and the jet then shuts down, this is an 
upper limit to the time taken by the cavity to move away from where it formed, 
hence to its age. In the notation above, = 2\Jr j g. As discussed, the early 
expansion of a cavity is likely to be supersonic, whereas motion of cavities under 
the control of buoyancy is almost invariably subsonic. Employing the simple 
compromise assumption, that a bubble is launched from the nucleus and travels 
at approximately the sound speed, it follows that the time it takes to rise to its 
projected position is then the sound crossing time, = R/cg. The sound speed is 
given by Cg = kT / {fimn) ~ llOOTg^^^ km s~^, where the ratio of specific heats, 
7 ~ 5/3, and fj, ~ 0.62. Generally, tc < ^buoy < U for observed cavities. 

Which timescale best approximates the true age depends on severa l unknown 
facto rs. Viscous stresses ( Revnolds et al. 20051 ) and magnetic stresses ( De Young 
20031 ) can reduce the terminal speed, so that the estimate given here may be 



high, causing ages to be underestimated and cavity powers to be overestimated. 
However, equating the age to the buoyant rise time presumes that cavity dynamics 
are controlled entirely by buoyancy. There are at least a few instances where this 
is clearly not so. 

Other unknown factors include the trajectory with respect to the plane of 
the sky and the effective location at launch. Projection causes rise times to be 
underestimated, but generally by less than a factor of two. If cavity dynamics 
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are controlled by the jet, all of the age estimates above are too long and powers 
would likewise be underestimated. Although adiabatic cavities expand as they 
move outward, the ratio r/R\s invariably a decreasing function of i?, so that the 
relatively large observed values of r/R (Figure [5]) add weight to the argument 
that most observed cavities have been caught early, before, or soon after becoming 
buoyancy dominated. Thus, outburst ages may generally be overestimated rather 
than underestimated. 

For example, consider the outburst in Hy dra A, for w hich shock models give 
a robust age estimate of ~ 1.4 x 10^ year ( Nulsen et al.. 2 005b'). The buoyant 
rise time for its large northern cavity (cavity E of Wise et al. II2007I ) is 5.2 x 10® 
year, some three to four times longer. This discrepancy would be expected if 
the dynamics of the cavities in Hydra A are still under the control of the jet. 
The same applies to most systems with shocks closely surrounding the cavities, 
and it is likely that shocks will be detected in deeper exposures of many more 
cavity systems. Because the shocks are relatively weak in most systems, their 
cavities cannot be advancing much faster than the speed of sound. Thus we get 
a crude estimate of cavity age from the sound crossing time. We note that this 
still overestimates the age of the northern cavity of Hydra A by a factor ~ 2. 



3.6 Stability of Cavities 

A significant issue is cavity disruption, which, in effect, can make cavities nonadi- 
abatic. Buoyant cavities are prone to Rayl eigh- Taylor and shear instabilities (e.g., 
Briiggen fc Kaiseil[200ll : ISoker et al.ll2002l ) , which can disrupt them, mixing all or 
part of their contents with the surrounding gas. The fate of a cavity's energy then 
depends on its constituents. Highly energetic particles can diffuse over large dis- 
tance s without depositing much of their energy as heat (e.g.. iBohringer &: Morfiii 
19881 ). Lower energy particles (such as nonrelativistic protons) will deposit energy 
locally. The fate of magnetic energy is unclear. Leakage of particles and fields 
may account for the radio mini-halo and extended nonthermal X-ray emission in 
the Perseus cluster ( Sanders et ahllioosl ). 



Abell 2052 is an example of a cavity system in the throes of disruption. The 
rim surrounding its n orthern cavity appe ars to be breaking apart where radio 
plasma is leaking out ( Blanton et al. 200 ll ). M87 is the nearest cavity system in 
a cluster ar id one of the best st udied. Its "bud" cavity is strongly suggestive of 
i nstability (Forman et al. 20051 ). However, the numerous small cavities in M87 
(|Forman et al.ll2007l ) suggest that cavities can fragment without being destroyed. 
Furthermore, the highly irregular shapes of the cavities in the simulations that 
demonstrate strong instability are at odds with the apparently regular shapes 
of observed cavities, also suggesting that real cavities are more stable than the 
simulated ones. 

Magnetic draping ( Lvutikov 20061 ) may help to explain this. A buoyant cavity 



can entrain magnetic field from the surrounding gas, stretching it and natu- 
rally creating a dynamically significant magnetic field around itself. Magneto- 
hydrodynamic simulations suggest that this effect, together with internal cav- 
ity fields, can help to s tabilize cavities ( De Youna 2003 : Robinson et al.l 2004 : 
Jones k de Younj l2005ll Growth rates of instab ilities are also affected by the 
viscosity ( Revnolds et al. 20051: Kaiser et al. l2005l ). which is poorly known (Sec- 
tion 1.5). As noted bv lPizzolato Sokei] (|2006l l. while a cavity is being inflated 
by a jet, its surface is often decelerating. Early in the lifetime of the cavity, the 
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deceleration can exceed the local acceleration due to gravity, preventing Rayleigh- 
Taylor instability. As discussed in Section 5.1, cavities that survive rising through 
several pressure scale heights, whether whole or as fragments, will liberate a sig- 
nificant part of their enthalpy as thermal energy in their wakes. 



3.7 Radio Lobe Composition Inferred from X-Ray Observations 
of Cavities 

Extragalactic radio sources are essentially bi polar outflows of magnetic field and 
relatiyistic particles ejected from an AGN (IBurbidgd Il956l: iBlandford &: Rees 
19741 : iBeeelman et al.lll984l : lOe Youndl200ll : iHarris k Krawczvnski 1120061 ^. Their 
structure includes a core associated with the AGN, oppositely collimated jets 
emanating from the core, and lobes that bloom from the jet terminals. Jets are 
narrow, collimated conduits that transmit mass, momentum, energy, and elec- 
tromagnetic field from the nucleus to the lobes, which in turn transmit much of 
the energy to the surrounding medium. Radio sources emit synchrotron radi- 
ation throughout the radio (and sometimes optical and X-ray) spectrum, from 
relativistic electrons gyrating along magnetic field lines. Synchrotron radiation 
reveals only the existence o f relativistic electrons and magnetic fields, not the ir 
momentum flux and power ( Begelman et al. 1984 : Harris &: Krawczvnski 20061 ). 
Charge neutralizing particles, such as protons for example, could carry most of 
the momentum without betraying their existence through the emission of radi- 
ation. We will focus on the lower synchrotron power FR I radio sources found 
in clusters, rather than the higher synchrotron power FR lis associated with 
powerful radio galaxies and quasars. Cavities and the associated shock fronts 
function essentially as gauges of the total energy output of jets, and they allow 
the contents of radio jets and lobes to be studied with greater certainty than 
was previously possible. These studies are revealing that even the jets of faint 
synchrotron sources can carry powers comparable to the luminosities of powerful 
quasars. 

The energy content of radio lobes is at least the sum of energy in particles, Ep, 
and magnetic field, Eb] 



Etot = EB + Ep = <^V^ + (1 + k)AL,ynB-^/\ 



(6) 



where V is the lobe volume, and $ is the volume filling factor of the magnetic field, 
B. The energy in electrons is related to the radio power per unit frequency, -Lgyn, 
as Ep, = ALsynB~^/'^ . The unknown factor, fc, accounts for unse en particle species 
()Dunn k Fabianll2004l : [Govoni k Feretti ll2004l : lDe Yound[200^ ) . so that the total 
particle energy can be expressed as .Ep = (l-l-fc).E'e. Corresponding approxim ately 
to the minimum energy condition ( Govoni k Feretti 20041 : De Youna 20061 ) . it is 
customary to assume energy equipartition between magnetic field and particles, 
yielding 



B, 



eq 



47r(l + k)AL 



syn 



2/7 



(7) 



For electron Lorentz factors of a few thousand, synchrotron radiation at radio 
frequencies requires ~ /iC magnetic field strengths. If the magnetic field strength 
greatly exceeds the equipartition value the magnetic field would control the jet 
dynamics. If the field strength is substantially below the equipartition value, the 
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particle pressure controls the dynamics. No models of magnetic field-dominated 
jets have been found to be stable over decades of scale, from the Schwarschild 
radius to the lobes, and thus it seems ui ilikely that the e nergy density in magnetic 
fields greatly exceeds that of particles ( De Youn^ 2006 1). 

The internal energy of the cavity, Etot, as well as its filling factor, and mag- 
netic field strength, B, can be constrained by X-ray observations. As discussed 
in Section 5.1, the internal energy ranges from pV for a magnetically dominated 
cavity to 3pV for a lobe dominated by relativistic particles (this is the enthalpy 
minus the work of inflation, pV). The gas pressure surrounding a cavity pro- 
vides a measure of the energy density in fleld and particles required to support it 
against collapse. Several studies have shown that at equipartition, lobe pressures 
are a pproximately an order of magnitude snia,ller than the sur rounding gas pres- 
sure ( Blanton et al. 2001 : Fabian et al. 2002 : De Youn^ 20061 ). One way around 
this is to suppose that the filling factor of the radio plasma, is less than unity. 
In that case, if the radio plasma is in equipartition at the external gas pressure, p, 
its filling factor would be $ = {peq/py^^, where peq is the equipartition pressure 
determined assuming a filling factor of unity. Thus, the filling factor of the radio 
plasma might be only percents and total lobe energies reduced by similar factors. 

However, this possibility has been ruled out by X-ray observations of cavities 
and shocks. X-ray count deficits over lobes are typically consistent with the X-ray 
emitting gas being completely excluded from radio lobes. While it is difficult to 
place stringent limits on the amout of X-ray emitting gas within lobes, cavities 
would not be evident in X-ray images unless the lobes displace most of the X-ray 
emitting gas. Furthermore, the expanding radio lobes are the pistons that drive 
shocks into the surrounding gas. To obtain significant shocks, the work done by 
an expanding lobe, / pdV, must be comparable to PsVs, where ps is the preshock 
pressure and Vg is the volume encompassed by the shock. For known shocks, this 
generally requires the lobes to displace most of the gas within the volume they 
occupy, i.e., it demands filling factors, ^, close to unity. The pressure support 
c ould then be supp lied by some combination of thermal or relativistic particles. 



De YoungI (|2006l ) modeled radio jets as pipe-like conduits of energy collimated 



by the surrounding gas pressure. The energy fiux was inferred from the pV 
work and buoyancy ages, based on X-ray observations alone, while the lengths 
and cross s ection s of the jets were taken from high-resolution radio observations. 
De Yound tood ) found that the energy in pV work alone is so large that the 



jets would decollimate unless most of the energy and momentum are carried by 
cold, heavy particles (e.g., protons) that do not contribute significantly to the 
internal isotropic pressure of the jet (but see Dunn, Fabian & Celotti (2006) for 
a discussion of electron-positron jets). De Young's results are consistent with the 
high ratio of relativistic particle energy to electron energy (A;) found by Dunn, 
Fabian, & Taylor (2005). Despite the many unknowns concerning, for example, 
jet stability, confinement, and acceleration mechanisms, De Young's analysis sug- 
gests that electrons are unable to supply the observed jet power alone and must 
be aided by heavy particles, or perhaps Poynting fiux. 

The fluid supporting th e cavities could be dominate d by a hot thermal plasma, 
or cosmic ray pressure (|Mathews fc Brighentil 120071 '). Constraints on such a 
fluid can be placed by asking what combination of gas temperature and den- 
sity would be required to provide pressure support without violating the X-ray 
surface brightness constraints. This technically challengin g measurement has 
yielded constraints for a few systems of kT > 15 - 20 keV ( Nulsen et al. 20021 : 
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Blanton et al]l2003l ) 



Sunyaev-Zeldovich (SZ) measurements in the submillimeter band in principle 
provide a novel means to discr iminate between a thermal and nonthermal cav- 
ity fluid (IPfrommer et alJl2nn5l V The SZ decrement is sensitive only to thermal 



gas. Thus the presence or absence of a decrement toward the cavities them- 
selves would provide a test. The effect is subtle and the observations are difficult 
to perform using existing instrumentation, but sh ould be possible in th e future 
with the Atacama Large Millimeter Array (ALMA: IPfrommer et al] 20051 ). These 



measurements coupled to deeper X-ray observations have the potential to place 
valuable new constraints on the composition of r adio jets. In the longe r term, 
SZ effect can also be used to detect shock fronts (I Cavaliere Lapi 20061 ) . Like- 



wise, GLAST may place interesting constraints on the cosmic ray content of 
caviti es through observations of t he pion decay continuum produced in cavity 
walls (jMathews Brighentil [20071 ) . 



3.8 Radiative Efficiency of Radio Sources 

The plot of jet power determined from X-ray cavity data versus total synchrotron 
power (core plus lobes) in Figure [7] shows a trend, as one might expect. The 



medi an ratio of jet (cavity) power to synchrotron power is ^ 100 (iBirzan et al. 



20041 ). However, the mean is much larger ~ 2800 owing to the large scatter. When 



considering radio flux from the lobes alone, the average ratio rises dramatically 
to ~ 4700 (L. Birzan, private communication). Figure [7] clearly demonstrates 
that most of the jet power is deposited into the surrounding medium and that a 
negligible fraction is radiated away. Although this result is not terribly surprising 
in a qualitative sense, we now have a quantitative measurement of this ratio 



that, interestingly, is on the high side of theoretical expectations (e.g.. iDe Young 
I2OOII ). Furthermore, the large scatter in the relationship, ranging between unity 
and several thousand, shows that synchrotron luminosity is a poor measure of 
true jet power. Factors contributing to the scatter include variations in age, 
field strength, and jet composition, but how these variables combine to create 
the scatter is not understood. Several radio faint cD galaxies with jet powers 
that equal or exceed the energy output of powerful quasars and radio sources like 
Cygnus A have been identified through X-ray observations, yet they would not 
have been identified as such using optical and radio observations alone. X-ray 
observations have revealed that under some circumstances black holes produce 
powerful mechanical outflows with little accompanying radiation. 

It is worth noting that adopting a pV-hased standard for AGN energy output 
blurs the canonical separation between high-power FR II radio sources and lower- 
power FR I radio sources that are typically found in the centers of clusters. 
Cygnus A is the best known FR II radi o source, and is th e most powerful radio 



source in the 3C cata log within z = 1 (cf. Young et al.l200j ). Seven objects in the 



Raffertv et al. ( 20061 ) sample of cavity systems exceed Cygnus A in mechanical 



power, yet none of them are considered to be powerful radio sources. The most 
distant object in the sample lies at z ~ 0.5, but most lie within z ~ 0.2. 

3.9 Simulations of Mixing by Radio Lobes 

The rising abundance gradients found in the cores of cD clusters are presumably 
established over a few gigayears by stellar evolution in cD galaxies. The gradients 



Heating Hot Atmospheres with AGN 



23 



should in principle then be sensitive to erasure by mixing induced by merger and 
AGN activity over that time frame, thus providing interesting constraints on both 
the outburst and merger history, as well as jet entrainment models. 

AGN-induced mixing has been e xplored recently by severa l group s using 2D 
and 3p hydrodynamical sim ulations (jBriigge lOmma et al. 200 4: lHeath et alJ 

200fil : iRoediger etla]l200fil ). The studies generally assumed an initial metallicity 
gradient added to a /3-like model atmosphere with gas particles tagged by metal- 
licity. Moderate jet powers of 1.4 x 10^^ erg s~^ to 6 x 10^^ erg s~^, and quasar- 
like jets exceed i ng 10 ^^ erg lasting 100 Myr to 3 Gyr have been explored. 
Roediger et aD (|2006l 'l tailored their simulation to the conditions in Perseus and 
included a prescription for metal injection from the central galaxy. 

These studies found that lower power jets allowed to run for 100 Myr or so 
produced only modest dredging that re duced existing met allicity gradients by 
at most a few tens of percent. However, Heath et al. ( 20061 ) found that quasar- 
like outbursts are able to reduce metallicity gradients to roughly 10% of their 
initial values by transporting the metals outward, primarily through convection 
and entrainment behind the cavities. Short duration (powerful) jets producing 
big cavities, and wide jets were most effective. Simulations generally produce 
anisotropic abundance distributions aligned along the jet prediction that 

can be tested using available X-ray and radio observations. However, these stud- 
ies considered jets launched at constant aspect angles into static atmospheres. 
An isotropic distribution of metals could be pre served t hrough turbulen t mix- 
ing and gas circulation (e.g., Mathews &: Brighenti 2003: Heinz et al.l 20061 ) or by 
launching bubbles on random traject ories, perhaps thr ough jet precession, as has 



been observed in some clusters fe.g.. lGitti et al.ll2006l ). 



3.10 Observations of Outflows and Mixing 

The outward mixing of metal-enriched gas has been invoked to explain why cen- 
tral metallicity peaks are broader than the l i ght profiles of cD g alaxies (e.g., 
David et al.ll200ll : IPe Grandi fc Molendil I2OOII : iRebusco et~aD I2OO5I ) . Moreover, 
there are several striking examples of plumes and shells of metal-enriched plasma 
and cold clouds that were apparently dragged outward by radio jets and lobes 
advancing into the ICM. A metal-enriched shell of gas was found nearly 9 kpc 
from the nucleus of NGC 1275 in the Perseus cluster ( Sanders et al. 20041 ) . and 
a plume of cool. X-ray emitting gas extends several tens of kpc along the radio 
axis of Hydra A (|Nulsen et al.ll2002l ;i. A metal-rich shell or cap of gas was found 
near the edge of the souther n cavitiy, 34 k pc from the nucleus of the central 
galaxy in the HCG 62 group (|Gu et alJl2007i ). However, the absence of a similar 
feature associated with the northern cavity and other circumstantial evidence for 
merger activity lead the researchers to propose the gas may have been stripped 
from an interloping galaxy. Arcs and s hells of Ha eraission surround the cavities 
and radio sources in several clusters (Blanton et al. 2001 : Fabian et al.l 2003bl : 



Crawford et al.l 20051 : Hatch et al. 2006 ), and cold mo lecular gas surrounds the 



inner cavities in the Perseus cluster ( Salome et al. 20061 ). The cooling time of the 



hot gas near the cavities is too long for the gas to have cooled locally, so it was 
probably dragged there from below. 
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4 STABILIZING COOLING FLOWS BY FEEDBACK 



4.1 Heating and Cooling Rates in Clusters 

The failure to f ind large quantities o f cooling gas with the expected properties of 
a cooling flow ( Peterson et alj 2003 ^1 implies that more than 90% of the energy 
radiated away is being replenished. Only a few percent of the gas associated with 
the cooling flow forms stars and even less accretes onto the central supermassive 
black hole. For an AGN to be a viable agent, it must be powerful, persistent, an 
efficient heater, and it must distribute the heat throughout the cooling region. If 
jets are underpowered, the remaining issues are moot. 

As discussed in Section 2, in terms of the classical cooling flow model, the 
power radiated from the core of a cluster can be expressed as Lx — 1.3 x 
W'^'^T^M2 erg s~^, where the classical cooling is given by M2 = M/(100 Mq year~^ 
This radiated power is equivalent to ~ 10% of the gravitational binding energy 
released by only ~ 0.02 Mq year~^ of accretion onto a black hole. The fact 
that AGN outbursts are so frequently associated with cool and cooling X-ray 
atmospheres suggests that the rates of cooling and AGN heating may be thermo- 
statically controlled. Three lines of evidence draw us to this conclusion. First, 
central cooling times are as short as ~ 3 x 10^ year, ~ 1/30 of t he Hubble time, 
in m any clusters and they are even shorter in elliptical galaxies (IVoigt &: Fabian 
iooi). The ~ 100 Mq year of cooling gas expected in the classical cooling flow 
model is rarely observed, implying that the energy radiated is being replenished 
on a sho rter timescale than t he small cooling times found within ~ 10 kpc of the 
nucleus ( Nulsen et al.ll2005bl ). 



If the heating rate were not tied to the radiative cooling rate, then it must 
exceed it to ensure that very little gas cools to low temperatures in the majority 
of cooling flows. This would only be possible if the heating rate exceeded typical 
cooling rates, giving rise to net heating of the ICM. Net heating would then drive 
up the central entropy, and in the long-term, the central cooling times would ap- 
proach the ages ~ Hq^ of most clusters, in contradiction of observations. Thus, 
it is difficult to maintain a significant proportion of clusters, groups, and galaxies 
with short central cooling times, unless heating rates are tied to cooling rates 
by some feedback loop. Second, the trend between cavity-based power estimates 
and the X-ray luminosity of cooling gas shown in Figure [HI strongly implicates 
AGN: they apparently know about the cooling gas and vice versa. Third, the 



implying entropy (energy) input at a level of 



(Donahue et al. 2006) 


cm^ ( 


McCarthv et al. 



2004l : rVoit &: Donahuell20()5l ). which is consistent with the observed level of AGN 
energy input. The flattening is seen in both radio active and radio quiet clusters. 
It suggests the operation of an intermittent heating mechanism that maintains a 
roughly steady power when averaged over tim es comparable to t he central cooling 
time. The coincid ent rates (70%) of cavities ( Dunn et al.l 20051 ) and radio emis- 
sion ( BurnsI 199d ) in cooling flow clusters further implic ates AGN, although heat 
conduction or other agents may a, ugment AGN heating (jRuszkowski &: Begelman 
2OO2I : iBrighenti fc Mathewsll2003l ). 

A quantitative comparison between X-ray (cooling) luminosity and jet power 
(heating) averaged over the lifetime of a cavi ty pW^huh is shown i n the heating 
versus cooling diagram (Figured]) taken from lRaffertv et al.l (l200fil ). In this dia- 
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gram, the cooling luminosity is the radiated power that must be replenished by 
heating. This quantity is found by subtracting from the total X-ray luminosity 
within the cooling radius the luminosity from gas that could be condensing out 
without violating observations. This correction amounts to less than 10% of the 
total X-ray luminosity. The diagonal lines in Figure [8] represent equality between 
heating and cooling assuming pV, ApV , and IQpV of heat input per cavity. Cav- 
ities filled with a nonrelativistic monoatomic gas ("hot bubbles") would supply 
approximately 2.5pV per cavity, whereas cavities filled with relativistic gas would 
supply roughly ApV per cavity (Section 5.1). The true effective energy per cav- 
ity could be substantially greater than 4:pV if much of the outburst energy has 
been dissipated by shocks (Section 5.4), or if other cavity systems exist below 
the threshold of detectability. Thus, the data points should be treated as lower 
limits. 

More than half of systems with detectable cavities liberate enough energy to 
balance or exceed radiation losses at the present time. Other systems including 
those without detectable cavities do not. Either this implies that AGN require 
help from other energy sources in order to suppress cooling, or it reflects the 
elusive nature of X-ray cavities and the transience and variable power output 
of AGN. If AGN outbursts are transient, then the time averaged AGN heating 
power needs to match the power radiated. Additional e nergy supplied by AGN 
associated with the broader cluster galaxy population ( Nusser et al.l 20061 ) will 
heat the gas beyond the cooli ng region but is apparently inconsequential within 



the cooling region itself (e.g., iBest et al.l 120071 ) 



The trend in Figure [8] does n ot take account of the 30% of cooling flow clusters 
that lack identifiable cavities ( Dunn &: Fabian! l2006l ) . which would in principle 
populate the lower right-hand side of the diagram. It is not known whether such 
objects have had AGN activity in the recent past, but the associated cavities and 
shock fronts are not visible in existing data, or whether other heating mechanisms 
are at wor k. Given that deepe r observations have invariably revealed more struc- 
ture (e.g., Forman et al. 200?! ). there is good reason to suspect that the fraction 
of clusters for which the AGN power is sufficient to balance radiative losses will 
grow beyond half as deeper Chandra images become available. 

Regardless of this, the existing data already s ugges t that AGN heating can 
balance cooling. For the sample of Rafferty et al. ( 20061 ). the mean cooling power 
is 6.45 X 10'^'' erg s~^, while the mean cavity power is 1.01 x 10'^^ erg (using 
ApV per cavity and the bouyant lifetimes). Assuming that cooling fiow clusters 
without evident cavities have similar cooling powers and zero cavity power, we 
can correct t he ratio of mean cavi ty power to cooling power for these with a 
factor of 0.7 (lOunn fc FabianI l2006l ^ . giving 1.1 for this ratio. Thus, within the 
substantial uncertainties, it is plausible that time averaged AGN heating powers 
balance radiative cooling. A more accurate assessment of the ratio of heating 
to cooling must await a study of a complete, unbiased sample of cooling fiow 
clusters. 

If we assume for the moment that all cooling fiows are suppressed by AGN, 
the scatter in Figure El may be a consequence of variable AGN power output. 



Nipoti &: Binnevi (120051 ) modeled AGN power output as a Gaussian process with 



a log-normal distribution at a fixed cooling luminosity, and an observationally 
motivated outburst timescale of 10* years. Their model implies that in any given 
system there is a good chance of finding smaller than average jet powers, because 
much of the power is generated by less frequent but more powerful outbursts. In 
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this context, the objects falhng below the 2pV - 4pV hnes may be in a lower than 
ave rage outburst state and thus may be in or moving toward a cooling cycle. In 
the Nipoti &: Binnev ( 20051 ) model, the powerful, rare outbursts are experienced 
by all systems. These outbursts may be responsible for the high entropy ped estals 
observed in radio quiet cD galaxies ( Donahue et al. 20061 ; Pope et al. 20061 ) . 

How effective AGN heating is over the lives of clusters depends on how and 
how much their power output varies over time. Although AGN power output 
is a strong function of halo gas mass, power outputs vary widely at a given 
mass. The systems with the largest cavities represent the extreme in power 
output. We do not know whether they are unique to some clusters or whether 
all systems occasionally experience them. These rare but powerful outbursts can 
easily dominate a cD galaxy's AGN power output over the age of a cluster. 

To summarize, under the assumptions outlined above, AGN are powerful enough 
to supress cooling in many and perhaps all cooling flow systems. However, this 
conclusion depends on how well the X-ray method traces true jet power and how 
efficiently cavity enthalpy and shock energy is converted to heat. 



4.2 Heating and Cooling in Elliptical Galaxies and Groups 

The cooling time of the hot gas in the centers of some gE galaxies is less than 
~ 10^ year, which is shorter than found in the cores of clusters. As in clusters, 
the hot gas there is expected to cool and form stars in the host gEs, yet it fails to 



do so at the expected rates of ~ 1 Mq year ^ ( Mathews &: Brighenti 20031 ). Like 



clusters, many gEs have cavities or other disturbances in the ho t gas implicatin, 
AGN as significant heati ng agents (jMathews Sz Brighentil l2003l : iDiehl Static: 
200?! ). Jones et al.l ( 2007 ) have assembled Chandra X-ray observations of a nearly 
complete sample of roughly 160 nearby gE galaxies, 109 of which show signifi- 
cant diffuse e i nissio n from hot gas. Of those, 27 have significant AGN cavities. 
Nulsen et al.l (2007) evaluated the relative rates of AGN heating and radiative 
cooling for the 27 ellipticals with detectable cavities, giving the results shown in 
Figure [9l Apart from minor differences, this figure can be viewed as extending 
Figure [8] to less massive halos. Although heating matches or exceeds cooling in 
roughly half of the cluster systems in Figure El AGN power exceeds radiative 
losses in all but one or two of the nearby gEs shown in Figure [H Roughly one 
quarter of Jones's gEs with significant emission from hot gas have detectable cav- 
ities, a detection rate that is similar to the overall rate of detection in clusters, 



but i s smaller than the 70% detection rate in cluster cooling flows (iDunn et al. 
iooi). Based on central cooling times, all 109 of the gEs with significant hot 



atmospheres are cooling flows. Assuming IpV per cavity, the total heating power 
of the 27 nearby gEs harboring cavity systems is ~ 2 x 10^^ erg s~^. This is to 
be compared to the total cooling power of ~ 10^^ erg s~^ for all 109 nearby gEs 
with significant emission from hot gas. Allowing for a cavity enthalpy of 4:pV for 
relativistic gas and the significant boost to this owing to shock energy, overall 
heating and cooling powers for the nearby gE sample match reasonably well. Al- 
though there is still significant systematic uncertainty, it seems that intermittent 
AGN outbursts are a plausible mechanism for preventing X-ray emitting gas from 
cooling and forming stars in nearby gE galaxie s. A simila r resul t was found indi- 
rectly for gEs culled from the Sloan survey bv lBest et aP (|2006l ). who determine 
jet power by assuming a one to one c orresp ondence with radio power, calibrated 
using cavity data from Birzan et al. ( 20041 ). 
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Because the binding energy per particle is lower in groups, the same nongrav- 
itational heating energy per particle will have a greater effect in groups than in 
clusters. Rosat studies of poor clusters and groups reported an "entropy floor" 
or pedestal in the gas entropy profile (e.g., Ponman et al.l 19991 ). It was argued 
that the entropy floor caused the steepening in the luminosity versus tempera- 
ture relation, which was more pronounced at group masses. Entropy floors can be 
produced by a number of processes, includin g AGN ou tflows, supernova driven 
galactic winds, and by the effects of cooling (iVoit 20051 ). With higher resolution 
data from Chandra the situation has become murkier. Entropy floors are no 
longer seen, but there is an apparent exces s of entropy at larger radii, beyond the 
inner cooling region (jPonman et al.ll2003l l. iMushotzkvl (|2004l l has argued that 
the steeper than self similar power law scaling of luminosity versus temperature 
for rich clusters extends to groups, but with a larger scatter at the group scale. 

Unfortunately, AGN heating is not as well studied in groups as in clusters. We 
know that the gas in many g roups is being distur bed by radio sources emanating 
fro ni their central galaxies (Croston et aP 120051 1. for example, as seen in HCG 
62 (|Morita et alll2006l : iNakazawa et al.ll2007l l. The cavity power of the lobes is 
sufficient to quench t he cooling flow but the over all significance of AGN heating 
in groups is unclear (iDwarakanath &: Nathll2006l l. 

The effects of AGN heating on intragroup gas has been explored by separating 
X-ray bright systems with central galaxies into radio-loud ari d radio -quiet bins. 
Analyzing an X-ray-selected sample of groups, ICroston et al.l (120051 ) found that 
in 63% of groups with a dominant central elliptical galaxy the dominant elliptical 
harbors an active radio source. Many of these radio sources are interacting with 
the hot gas filling the groups. ICroston et all (|200,5l l found that the radio loud 
groups are slightly hotter on average at a given X-ray luminosity than radio quiet 
groups, which they attribute to AGN heating. They point out that AGN heating 
at the level they find extending over 5 x 10^ year would be able to supply the ~ 1 
keV p er particle of excess energy required to preheat groups (e.g. Voit &: Donahuel 



20051) 



However, in a Chandra study of 15 groups, I.Tetha et al.l (|2007l l found only 
modest steepening of the gas temperature in radio-loud groups compared to radio- 
quiet ones, and steeper entropy gradients in groups with brighter central galaxies 
and presumably more massive nuclear black holes. It is unclear whether these 
trends implicate AGN heating or heating by the galaxies themselves. Jetha et 
al. found power law entropy profiles extending to small radii, and only small 
differences in the gas profiles of radio-loud and radio-quiet groups. They found 
that AGN heating may have a more significant effect in quenching cooling than 
for preheating. 



4.3 Uncertainties in Estimating Jet Power 

Model assumptions about the time and energy dependence of outbursts are a 
significant source of uncertainty in the calculation of mean cavity (jet) power. If 
AGN inject energy in a series of short, isolated outbursts, the interval between 
outbursts characterizes tbub- This interval is impossible to estimate for individual 
systems in the absen ce of multiple generations of cavities, such as those in Perseus 
( Fabian et al. 200d ) and Abell 2597. Multiple cavity systems are also required 
to estimate the mean energy per outburst. The ghost cavities in Perseus and 
Abell 2597 imply an outburst every 60 - 100 Myr or so, and in both clusters 
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the earher outbursts were stronger than the current ones, indicating variabihty 
in outburst strength. The ripples in Perseus indicate that some outbursts may 
occur on shorter timescales than the cavity ages, implying the outburst period is 
also variable. 

If outbursts are more nearly continuous, i.e., the jets remain active for a sub- 
stantial fraction of the outburst cycle, then the current jet power is a more useful 
measure of the average value. In general, the off-time must be included to measure 
a population average. The large systems ( Hydra A, MS 0735. 6+7421, Hercules A) 
appear to be operating in this mode (e.g.. Wise et al. | [2007). Most observational 
treatments implicitly assume a more nearly continuous distribution of outbursts. 

Sources of uncertainty in the measurements of cavity energy include the volume 
estimates from the projected cavity sizes and shapes and the unknown composi- 
tion of the cavity plasma, which combined give an uncertainty of at least a factor 
of several. In addition, the energy per outburst is probably underestimated ow- 
ing to adiabatic los ses, cavity disruption, undetected cavit ies, and the omission of 
shock energy (e.g., iNusser et aDl2006l : iBinnev et aT]l2007l l. These effects are off- 
set to some degree by the unknown fraction of outburst energy that is converted 
to heat. 



5 HEATING MECHANISMS 
5.1 Cavity Heating 

The energy required to create X-ray cavities around radio lobes is the sum of the 
pV work required to displace the X-ray emitting gas and the thermal energy of 
the contents of the lobe, i.e., the enthalpy, 

H = E + pV = ^^pV, 

where p is the pressure in the lobe and V is its volume. The second form applies 
if the lobe is filled with an ideal gas with constant ratio of specific heats, F. If 
the lobe is dominated by relativistic particles, F = 4/3 and H = 4:pV, whereas if 
it is dominated by nonrelativistic gas, T = 5/3 and H = 2.5pV. Lobes may also 
be dominated by magnetic field, in which case H = 2pV. Other possibilities lie 
between these extremes, so that, although the equation of state (F) for lobes is 
not known, lobe enthalpy is likely to fall in the range 2pV - ApV. 

Simulations show that buoyant cavities can heat the surrounding gas as the y 
rise through a cluster atmosphere (jBriiggen fc Kaiserl[2nni : lEevnolds et"aDl2nn2l ^ . 



Churazov et al. I (120021 ^ argued that essentially all of the enthalpy of a rising cavity 



can be thernialized in its wake and a simple mechanism for this is outlined in 
( Birzan et al. 20041 ) . As a buoyant cavity rises, some X-ray emitting gas must 



move inward to fill the space it vacates (Figure flOl) . so that gravitational potential 
energy is turned into kinetic energy in the ICM. The potential energy released as 
a cavity rises a distance 6R is 

5U = Mg5R = Vpg5R = -V^5R = -V6p, (8) 

where M = pV is the mass of gas displaced by the cavity, V is its volume, p is the 
density of the surrounding gas and g is the acceleration due to gravity. The third 
equality relies on the surrounding gas being close to hydrostatic equilibrium. 
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so that pg = —dp/dR. The last equahty expresses the result in terms of the 
change in pressure of the surrounding gas over the distance 5R^ (dp/ dR) 6R = 6p. 
Because the rising cavity moves subsonically ( Churazov et al. 200 ll ). its pressure 
remains close to that of its surroundings and we may regard this as the change 
in its pressure. The first law of thermodynamics, dE = TdS — pdV, expressed 
in terms of the enthalpy is dH = TdS + Vdp. Entropy remains constant for an 
adiabatic cavity (radiative losses from the cavity are negligible) , so that this gives 
6H = V 6p. Thus, Equation [8] shows that the kinetic energy created in the wake 
of the rising cavity is equal to the enthalpy lost by the cavity as it rises. 

Regardless of the viscosity, we should expect this kinetic energy to be dissi- 
pated, creating heat locally in the wake of the cavity. If the viscosity is high, 
i.e., the Reynolds number is low, the motion is damped viscously in a laminar 
wake that is comparable in size to the cavity. If the Reynolds number is high, 
the wake is turbulent and kinetic energy is damped on the turn-over timescale 
of the largest eddies (also comparable to the size of the cavity). The turbulent 
cascade maintains the dissipation rate by propagating energy to sufficiently small 
scales for viscous dissipation to keep pace (in fact, the turbule nt spectrum is fully 
characterized by the dissipation rate; Landau &: Lifshi"t3 198?! ). In both cases, ki- 
netic energy is damped before diffusing far from the axis on which the cavity 
rises. This leads to the important result that the enthalpy lost by a buoyantly 
rising cavity is thermalized locally in its wake, almost regardless of the physical 
properties of the cavity and the surrounding gas. 

Heating by cavity entha l py is the physical basis of th e "effervescent" heat- 
ing m odel ( Begelman 2001 : Ruszkowski &: BegelmaiJ 20021 : Roychowdhury et al 
20041 ). Chandra observations have established not only that cavities are common 



but also that many clu sters contain multip le cavities, though not all coincident 



with active radio lob es (Fabian et al. l200d 'l. The numerous cavities seen in deep 



observations of M87 ( Forman et al.l 200?! ) also fit this model. The enthalpy of an 
adiabatic cavity depends on the surrounding pressure as H\y = H]^fi{p/pQ)^^~^^^^ , 
where H]^^ is the initial enthalpy of the cavity and po is its initial pressure. If the 
mean power injected by an AGN as cavity enthalpy is Lb, then the mean heating 
rate per unit volume averaged over the sphere at radius R due to liberated cavity 
enthalpy is 



Lh d 
' AttK^ dR 



p\ 
Po) 



{r-i)/r 



(9) 



where pq is the pressure at the radius, i?o, where the cavity is formed. The 
value of i?o depends on outburst details, but it is always nonzero. Very little 
cavity enthalpy is thermalized wit hin Rn- W i th mi nor modifications, this is the 
efferves cent heating rate used b y Begelman ( 200ll ) and successors. Using this 
model, Ruszkowski fc Begelman! ( 2002 ) argued that cavity enthalpy needs to be 
supplemented by thermal conduction from the surrounding cluster in order to 
prevent cooling. However, iRoychowdhury et al. (2004 ) find that this conclusion 
depends on the assumed total rate of AGN heating. Note that these are all 
ID models, employing mixing length theory to model convection. Such models 
probably omit significant features of the full, 3D flow. 

As well as providing better models for the convection associated with anisotropic 
cavity heatin g, 3D simulations al low more realistic tr e atmen ts of cavity stability 
and mixing dOuilis et al.]l200lh . iDalla Vecchia et al. (j2004] ) made 3D hydro dy- 
namic simulations for the longer-term effects of cavity heating by injecting energy 
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to simulate the formation of radio lobes at random positions with a Gaussian 
distribution around cluster centers. Although they found that the resulting heat 
input can prevent catastrophic cooling, their model clusters do not produce cool 
cores as observed. They argue that additional preheating, taking place earlier 
in the collapse hierarchy, could resolve this disagreer nent. AGN heating rates in 
this model are fixed, not determined by feedback; and IPalla Vecchia et ^il 
conclude that prevention of catastrophic cooling is not sensitive to AGN power. 
A significant channel of heating for their model is the mixing of gas from the 
simulated cavities with the general ICM. Relativistic particles can be slow to 
transfer energy to gas (e.g., Bohringer &: Morfilll 19881 ) and the fate of magnetic 
fields is poorly understood, making it unclear how effectively radio plasma mixes 
with and heats the ICM in practice. 



5.2 Heating by Weak Shocks 

Although a substantial energy, comparable to the cavity enthalpy, is required to 
drive the weak shocks seen in association w ith AGN outbursts in some clusters 
( Forman et al. 20051 : McNamara et al. 20051 ). much of the shock energy ends up 
as additional potential energy in the gas. This helps to delay cooling by reducing 
gas density and increasing the total energy that must be dissipated. However, 
the key requirement on the process that prevents the gas from cooling is to make 
up for entropy lost by radiation from the gas. Shock heating probably p lays 



a significant role in this, especially close to the AGN ( Fabian et al. 20051 ). A 



fundamental distinction between sound waves and shocks is the entropy created 
by dissipation at shock fronts. For weak shocks, the entropy j ump per unit mass, 
Ag, i s proportional to the cube of the shock strength (e.g.. Landau &: Lifshita 



19871 : IPavid et al.ll200ll : measured here by the fractional pressure increase, 5p/p, 



where p is the preshock pressure and 6p is the pressure increase across the shock) . 
The equivalent heat input per unit mass is T AS, where T is the gas temperature. 
To lowest nonzero order, this equivalent heat input amounts to a mean heating 
rate per unit volume due to repeated weak shocks of 



(7 +l)ujp / dp 
1272 2^ 



(10) 



where 7 is the ratio of specific heats for the ICM and the interval between out- 
bursts is expressed as 2it/uj (compare the sound heating rate. Equation II 111. 



Th e repeated weak shocks around M87 seen by Forman and colleagues fiForman et al 
20071 ^ demonstrate the possible significance of weak shock heating. Based on its 
X-ray surface brightness profile, the innermost shock, 0.8 arcmin (4 kpc) from 
the AGN, has a Mach number of 1.4, so that its equivalent heat input is only 2% 
of the gas thermal energy. However, there is another comparable shock at about 
twice the radius, suggesting repeated outbursts every ~ 2.5 x 10^ year, while the 
cooling time of the gas is ~ 2.5 x 108 year. Over the cooling time, this leaves 
amp le time for m ultiple weak shocks to make up for the energy radiated by the 
gas ( Nulsen et al. 1[2007). In M87, weak shocks are capable of preventing gas near 
the AGN from cooling. 

When a shock overruns a cavity, the high sound speed of the cavity causes 
the shock to propagate faster through the cavity than around it. This results 
i n the formation of a vortex ring around the c avity after it has b e en ov errun 
(jChurazov et al.ll200ll : lEnfilin fc Briiggenl l2002l ^ . iHeinz fc Churazovl (|2005l ) note 
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that this process, an example of Richtmyer-Meshkov instability, can increase the 
fraction of shock energy that is thermalized in the ICM, especially for weak 
shocks, because it converts shock energy into localized kinetic energy that can 
then be dissipated as heat. If the ICM contains many small cavities, Heinz & 
Churazov find that the attenuation length due to the Richtmyer-Meshkov insta- 
bility is inversely proportional to the fraction of the ICM occupied by the cavities. 
For example, if cavities in the Perseus cluster have a filling factor of ~ 10%, most 
of the energy of the weak shocks would be dissipated within the central 100 kpc. 
The significance of this process is hard to assess, because small bubbles are not 
generally accessible to observation. If, for example, bubbles occupy ~ 10% of the 
volume of the ICM, then they also contain ~ 10% of its t hermal energy. Becaus e 
the bubbles rise at their terminal speeds, Vt ^ VK\/r/R ( Churazov et ahlboOll ). 
they must be replaced continually. The power required to maintain a large bubble 
filling factor is substantial, unless the bubbles are very small, and a large part of 
that power would be dissipated as bubble heating (Section 5.1). 



5.3 Heating by Sound Damping 



Repeated we a k shock s may also be regarded as a superposition of sound waves. 
Fabian et all (|2003al ^ showed that viscous damping of sound waves generated 
by repeated AGN outbursts may represent a significant source of heating. The 
he ating power per unit volume due to dissipation of a sound wave can be expressed 
as (iLandau k Lifshit j[T987l : iFabian et al.ll2005l ^ 



2u (7-i)^^r 

3 27P 



(11) 



where p, T, p and 7 are the density, temperature, pressure, and ratio of specific 
heats of the gas, respectively, k is the thermal conductivity, u is the kinematic 
viscosity {v = p/p, where p is the viscosity), uj is the angular frequency, and 6p is 
the pressure amplitude of the sound wave. This expression includes both viscous 
and conductive dissipation. Both terms in the leading c oefficient have the form 
of a mean free path times a thermal speed ( Spitzer 19621 ). For an unmagnetized 
plasma, the mean free paths of the electrons and protons are the same. For 
the kinematic viscosity, the thermal speed is that of the protons, whereas in the 
conductive term it is that of electrons. Thus conductive dissipation would be 
greater in the absence of a magnetic field. 

Heating by weak shocks is a separate mechanism from heating by sound dis- 
sipation. In nurnerical sir nulations, shocks ar e controlled by viscous stresses 
(|Ruszkowski et al.ll2004al lbl: iBriiggen et al.ll2005l '). so that these two processes are 
lumped together. Dissipation of sound depends on the transport coefficients, but 
shock heating does not. Therefore, uncertainty in the transport coefficients (Sec- 
tion 1.5) translates directly into uncertainty in the heating rate due to sound 
dissipation. If the transport coefficients are suppressed by no more than an order 
of magnitude, about the range sugge sted by observations, ana ly tical estimates 
( Fabian et al. 2003al ) and simulations ( Ruszkowski et aP 2004a bl: Briiggen et al 



20051 ^ agree that sound dissipation plays a significant role in converting AGN 
energy into heat in the ICM. The two mechanisms have distinct dependencies on 
the parameters (Equations [10] and [TT]) , leading to distinct spatial distributions of 
heating. 
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Only the fundamental frequency should be used in applying the sound dissi- 
pation rate (Equation [TT]) . The thickness of a weak shock front can be expressed 
as li; ~ A/(M — 1), where A is the effective particle m ean free path and M is 
the shock Mach number (e.g., Landau &: Lifshi"t3 19871 ). When the width of the 
shock front becomes comparable to or larger than the wavelength, the disturbance 
transitions from a shock wave to a sound wave. At this transition the sound dis- 
sipation rate, which depends on transport coefficients, and the shock dissipation 
rate, which does not, cross over. When the thickness of a shock is considerably 
smaller than the distance between successive shocks, the shock front contributes 
substantial power in high harmonics when decomposed into a Fourier series. The 
fundamental frequency is determined by the shock repetition rate. Because the 
dissipation rate varies as io'^ (Equation [TT]) the linear treatment suggests that 
the high frequency components will dissipate quickly, producing significant heat- 
ing. However, the entropy created at the shock front is determined entirely by 
the Rankine-Hugoniot jump conditions, regardless of the transport coefficients. 
Nonlinear effects adjust shock thickness to make the dissipation rate match that 
required by the jump conditions. Although shock thickness is determined by the 
transport coefficients, the dissipation rate for the high frequency terms is not. 
This dissipation is accounted for as the shock heating rate (Equation [T0|) . Only 
the dissipation rate of the component at the fundamental frequency depends di- 
rectly on the transport coefficients and only it should be regarded as subject to 
sound damping. 

In the absence of magnetic fields, the kinematic viscosity, v, and the coefficient 
of the conductive dissipation rate, nT/p in Equation [TT] both scale approximately 
with temperature as T^/"^. If this scaling applies in reality, the sound dissipation 
rate is sensitive to the ICM temperature, which would reduce its significance 
in cool systems. However, because the sound dissipation rate is also sensitive 
to frequency, higher outburst rates in cool systems could offset the reduction in 
transport coefficients. iFahian et al.l ^ ) find that outbursts occur every ~ 10 
year in the Perseus cluster, where the ICM temperature is ~ 4 keV, whereas the 
period of outbursts in M87 is ~ 2.5 x 10^ year for an ICM temperature 1—2 keV 
(iNulsen et alll2nn7l ). This hints that the feedback process may be able to adjust 



the frequency as well as the power of AGN outbursts. 



5.4 Cavity Enthalpy versus Shock Energy 

The division of jet energy between cavity enthalpy and shocks is affected signifi- 
cantly by the history of an outburst. Here, shock energy is used loosely to mean 
the total work done on the ICM by an inflating cavity, most of which ends up 
as thermal and gravitational potential energy in the hot gas. In fact, because 
cavity expansion lifts gas outward, it invariably reduces the gas pressure after a 
shock has passed. As a result, the gas expands and its total thermal energy is 
generally decreased after the outburst, especially when shocks are weak. Thus, 
the ultimate repository for most of the shock energy is probably potential energy 
in the ICM. 

At one extreme, consider a jet that dumps its energy explosively into the gas in 
a single brief event. Initially this would create a tiny cavity with volume, Vi, and 
pressure, pi, much greater than the pressure of the surrounding gas. The thermal 
energy of the cavity, Ei = PiVi/ (T — 1), would equal the energy deposited by the 
jet (radiative losses are assumed negligible throughout this process). This cavity 
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would expand explosively, driving a strong shock into the ICM, until its pressure 
reached that of its surroundings, pf <C Pi- Because the expansion is adiabatic, 
the final energy of the cavity is 

1 1 /«f\(^~^^/^ 

= = =^,(g) . (12) 

which can be much smaller than its initial thermal energy. The balance of the 
initial energy is the work done by the expanding cavity on the ICM, i.e., the shock 
energy. This extreme case illustrates that, in principle, there is no upper limit on 
the fraction of the energy from an AGN outburst that ends up in the ICM. At 
present, the known cavities in clusters are all weakly over pressured. However, the 
southwestern lobe of Cen A has a pressure th at is two orders of magnitude greater 
than that of the surrounding, unshocked gas (iKraft et al.ll2003l l. showing that the 



explosive extreme can be approached for systems in poorer environments. 

At the opposite extreme, consider a cavity that is inflated gently. Its pressure 
would remain close to that of the surrounding ICM throughout the expansion. 
To the extent that the pressure also remains constant during cavity inflation, the 
ratio of the work done by the expanding cavity to its final thermal energy would 
then be r — 1. In practice, the pressure will generally decline as a cavity inflates 
and expands outward into lower pressure gas. This causes additional adiabatic 
energy loss, boosting the ratio of the work done to its flnal thermal energy. 

The main conclusion here is that work done by an expanding cavity is rarely 
less than its thermal energy and may be considerably greater. This is broadly con- 
sistent with findings for observed clusters, t hat the cavity enthalpy i s comparable 



to the energy required to drive the shocks (iMcNamara et al.l l2005l : iNulsen et al 



This argument also shows that the relative energies contain information 



about the history of outbursts. 

5.5 Distribution of AGN Heating within a Cluster 

Ignoring all dissipation, we can use conservation of energy to estimate the radial 
dependence of the amplitude of sound waves and weak shocks, 

^OCp-l/2r-3/4i?-l, (13) 
p 

where p{R) and T{R) are the density and temperature of the ICM, and R is the 
radius. With this scaling, the ratio of weak shock (Equation 1101) to sound heatin g 



rate (Equation II ip almost invariably decreases with radius ( Fabian et al. 20051 ). 
As might be expected, this arises because the shock heating rate is more sensitive 
to 5p/p, which decrease s with the radius. The radial temperature gradient that 



typifies cool cores (e.g., lAllen et al.ll200lh probably decreases the ratio further. 



as the transport coefficients increase with radius. 

Because Equation [13] assumes that energy is conserved, it only applies after 
the expanding radio lobes stop driving the pressure disturbance, i.e., after the 
disturbance separates from the cavities. Using it to extrapolate to small R would 
overestimate the relative strength of the disturbance there. The assumption of 
spherical symmetry, which is adequate at large distances from the radio lobes, 
also fails at radii smaller than the lobe separation, where the off-center energy 
deposition must be taken into account. 
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The ratio of cavity heating (Equation [9]) to shock heating rate (Equation llOp is 
proportional to t!2.^V2-i/r2i5/4-i/r ^ p3/2-i/r2i9/4-i/r^ where the Kepler speed 

is vk = V'gR and g is the acceleration due to gravity, and the ICM is assumed 
to be hydrostatic. The second form relies on the further approximation, v'^ oc T. 
The relatively gentle temperature rise with radius in cool cores is insufficient 
to offset the density decrease, so that this is generally a decreasing function of 
the radius (for reasonable values of T). Thus, cavity heating is more centrally 
concentrated than weak shock heating. As noted above (Section 5.1), although 
it is formally the most centrally concentrated of these three heating processes, 
cavity heating is ineffective inside the radius where the cavities are formed. Thus, 
weak shock heating is likely to be the most significant heating process closest to 
the AGN. Cavity heating probably takes over this role irnmedi ately outside the 
region where the radio lobes are formed (|Voit k Donahue! l2005l l. 



If AGN outbursts deposit comparable amounts of energy in shocks and cavity 
enthalpy, then the reasoning of this section suggests that the dominant mode 
of AGN heating changes with radius. Closest to the AGN, weak shocks (or, 
possibly, sound dissipation) are likely to be most significant. Note that weak 
shock heating can plausibly stop the innermost gas from cooling in M87 (Se ction 
5.2). The total rate of shock heating may not be large ( Fabian et al. 20051 ) . but 



because the gas closest to the nuclear black hole is the most likely to be accreted, 
the heating process at work on that gas plays a critical role in any AGN feedback 
cycle. Cavity heating may well take over beyond the radius where the radio lobes 
are formed. On larger scales, sound damping may become the dominant AGN 
heating process. On even larger scales, thermal conduction can play the dominant 
role in the hotter clusters. In short, it appears likely that no single AGN heating 
process is the most significant. It may also be that AGN heating does not act 
alone to prevent copious gas deposition and star formation. 

5.6 Energy Injection by Radio Jets 

In order to understand the process of AGN feedback, it is necessary to understand 
how AGN outbursts are fueled and triggered (Section 7.3), a s well as the spa- 



tial d istribution and form of the energy deposited by jets (e.g., lOmma &: Binney 



20041 ) . Observations of shocks and cavities created by AGN outbursts have moti- 



vated the developme nt of increasingly sophi s ticate d simulations of the interaction 
of jets with the ICM. IVernaleo fc E'evnoldsl (|200fil ^ found that even a variable hy- 



drodynamic jet flowing into a static atmosphere is incapable of tranferring a 
significant fraction of its energy to the atmopshere. In order to avoid channeling 
jet energy beyond the halo, s ome additioii a l phy sics is required. 

A recent jet simulation by Heinz et al. ( 20061 ) is situated in a cluster atmo- 



sphere drawn from a realistic simulation of heirarchical structure formation. The 
simulated atmosphere includes the ongoing effects of cosmological evolution, such 
as merger shocks and turbulence. The jet also wobbles, simulatin g a dentist dri ll 
effect that is thought to result from instabilities along the jet (Scheuer 19741 ). 



Rather than punching through the atmosphere, the jet deposits much of its en- 
ergy into gas close to the AGN. This produces a fairly realistic simulation of 
Cygnus A (Figure [TT]) . Simulations of this sophistication need to be coupled 
with a realistic model for fueling and triggering outbursts in order to fully test 
feedback models. 

Although radio jets appear to be the main route by which AGN energy is 
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conveyed to the hot gas in nearby galaxies and clusters, AGN outbursts produce 
uncollimated winds and intense radiation that can also heat the gas. Powerful, 
uncollimated nuclear winds should couple strongly to the surrounding gas, making 
them an effective means of heating. However, they do not appear to be a majo r 



route for AGN energy injection in nearby systems (e.g., Krongold et al. 200?! ). 



Photoionization and inverse Compton heating can also c ouple a small fraction of 
the energy radiated in an AGN outburst to the hot ISM ( Ciotti &: Ostriker 200ll : 
bazonov et ai.ll2nn.^ ;i. Among other things, the significance of this process depends 
on the relative values of the g as temperature and the Compton temperature of 
the AGN radiation (~ 2 keV; ISazonov et all 120051 ') . making it most effective in 
lower temperature systems during the quasar era. Again, there is no evidence of 
such heating in nearby systems. 

6 HEATING WITHOUT FEEDBACK 



6.1 Conduction 

Whether or not thermal conduction is an effective heating agent has been a con- 
troversial issue for decades. Plenty of thermal energy is available in the outskirts 
of clusters to heat the cores. The controversy concerns whether and how it can be 
effectively tapped. Assuming the magnetic fields threading clusters are tangled 
(Section 1.4), the conductivity is suppressed below the classical Braginskii value, 
inhibiting inward heat flux (Section 1.5). The magnetic field topology is not ob- 
servable, so thermal conduction is poorly understood. Apart from the difficulties 
of maintaining local equilibrium betwee i i conductive heating and radiation (the 



fine-t uning problem; iStewart et al.lll984l : iMalagoli et al.lll987l : iBregman fc David 
19881 ). there are several clusters where heat conduction cannot balance coolin, 



throughout the core, even at the full B raginskii rate (jZakamska &: Naravanll200^ 
Voigt fc Fab"i^l2004l : IWise et al.1 [2003 1. 

On the other hand, if the thermal conductivity is only mildly suppressed (Sec- 
tion 1.5) , conduction may be effective in the outer reaches of the cooling regions of 
clust e rs (iNara van fc Medvedev 2001 ; Zakamska &: NarayanI 20031 : Voigt Sz Fabian 
20041 : IWise et al. 20o3)7 which can substantially reduce the power required from 
AGN heating. The heat flux depends strongly on temperature, as ~ 7"5/2^ qq 
is less important in l ower temperature systems, including gE galaxies. 

Dolag et all (|2004l l investigated the effect of including thermal conduction in 



numerical simulations of galaxy cluster formation. With the thermal conductivity 
set to one third of the Braginskii value, they found that although the temperature 
structure of the ICM is significantly modified, the fraction of baryons that cools 
and turns into stars is little affected. This supports the argument that thermal 
conduction acting alone does not explain why star formation is inhibited in cluster 
cores. 



6.2 Other Heating Mechanisms 

Many mechanisms have been proposed to prevent gas from cooling to low tem- 
perat ures at the c enters of cooling flows, but very few of these involve feedback 
(e.g., Fabian 19941 ). As discussed in Section 4.1, without feedback it is difficult 
to account for the many clusters with short central cooling times. Nevertheless, 
such mechanisms could significantly reduce the demands on AGN heating. Most 
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other mechanisms that have been proposed to stop the gas from coohng rely 
on the energy availa ble from mergers or, more generally, cluster binding energy 
^cf. lBlandfordll2003l l. Preventing the core gas from cooling in a system like the 
Perseus cluster requires ~ lO^'^ erg over a Hubble time. Because a major merger 
releases ~ 10^^ ergs, only a small part of that energy needs to be tapped in order 
to prevent gas from cooling and forming stars. Set against this, cool cores oc- 
cupy a small fraction, ~ 0.001, of the total cluster volume, presenting a relatively 
small target for undirected heating. Furthermore, the stably stratified cluster 
atmosphere, which resists the inward propagation of turbulence, combined with 
high central pressure (density), tends to shield the central gas from disturbances 
in the more tenuous gas that inhabits the bulk of a cluster. 

Simulations show that major mergers can disrupt some cluster cooling flows 
although, even for head-on mergers, the effect need not be long lasting (e.g. 



Gomez et al1l2002l l. Higher merger rates may explain t he reduced fraction o f 



bright clusters with strong central X-ray peaks at z > 0.5 ( Vikhlinin et al.ll2007l ) 
The processes that make cooling flow clusters different f rom clusters with Ion; 



central cooliiig tini es are still the subject of debate (e.g., iMcCarthv et al.ll200 



O'Hara et all l2006l l. but it seems likely that merger history plays some role. 



Thermal conduction should make it very difficult to re-establish cooling cores in 
hot clusters, unless it is suppressed by a large factor (Section 1.5). 

Simulations of the growth of structure consistently find that t he ICM should 
be tu rbulent, with typical turbulent velocities ~ 100 km s~^ (e.g.. lKravtsov et al 
Dissipation of turbulence could be a significant source of heating. Out- 



flows from central AGN may create a similar level of turb ulence, as suggested by 
observations of the Perseus cluste r (IChurazov et al.l 120041 . so that the source of 
turbulence in cool cores is unclear, bennis X. Cha.drani m used a semiempir- 
ical model to argue that a combination of thermal conduction, turbulent dissipa- 
tion, and turbulent diffusion could balance radiative losses for turbulent velocities 
in the range 100 - 300 km s~^. However, their model does not explain why the 
turbulence should have the spatial distribution or power required to balance ra- 
diative losses locally. 



7 FEEDBACK AND GALAXY FORMATION 



7.1 Star Formation in cD Galaxies 

Brightest clusters galaxies (which we refer to for convenience as cD galaxies) , with 
masses upward of 10^^ Mq and with halos extending hundreds of kiloparsecs 
into the s urround ing cluster, are the largest and most luminous galaxies in the 
Universe ( Sarazinl flOSS). cDs are similar in appearance to gE galaxies, but there 
are significant differences. Their central surface brightnesses tend to be lower, 
their velocity dispersions rise less steeply with increasing luminosity, and they of- 
ten possess stellar envelopes that lie above the R^^^ law profiles that characterize 
gE galaxies. Their locations at the centers of clusters suggest they grew to such 
enormous sizes by swal lowing stars and gas from i ieighboring ga laxies through 
mergers and stripping (jGallagher &: Ostrikeii Il972l : iMerrittI Il985l). This process 
is augjmented at late times b y the accretion of intracluster gas (jCowie &: Binnev 
19771 : IPabian k Nulsen|[T977l l. 

The stellar populations of most cDs are metal rich and dormant, and they 
rarely show optically active nuclei. Only 10% of optically selected cDs in non- 



Heating Hot Atmospheres with AGN 



37 



cooling flow clusters harbor nuclear ernission line fluxes abo ve an equivalent width 
of a few angstroms ( Best et aDbood : lEdwards et al. 2007I') . In comparison, 20% 



of galaxies lying outside of cluster cores (jBest et al. 



20061 ) have detectable line 



emission, suggesting that the gas processes driving emission line activity in most 
galaxies are suppressed in cluster cores. The detection rate and strength of neb- 
ular emi ssion in cD galaxies rises dramatically, to at least 45%, in cooling flow 
clusters (jCrawford et al.lll999l : [Edwards et al.ll2007l ). These emission line systems 
often extend tens of kiloparsecs from the nucleus ( Heckman et al. 19891 ). They 
appear to be powered by a combination of shock heating, stellar photoioniza- 
tion, and irradiation by the surrounding X-ray gas, but rarely by photoionization 
from an AGN. At least half of these systems are expe riencing star formation (e.g., 
Johnstone et al.lll987l : iMcNamara k 0'Connelll989l ') , and the nebular clouds and 
sites of star formation a re usually e mbedded in pools of 10^ - 10^^ Mq of atomic 
and molecular gas fe.g.. lEdgel[200ll ). The ratio of young stars to galaxy mass is 
typically only ~ 10~^. The corresponding star formation rates are typically a few 
M f;^ year~^, but in extreme cases they appro ach or exceed 100 year~^ (e.g., 



Crawford et al.l 1999 : McNamara et al. 2006l l. Young stellar population masses 
of 10® Mq to 10^" M0 are found in these systems. The largest starbursts ri- 
val those observed during the most rapid period of galaxy growth at z = 2 - 3 
( Juneau et al. 20051 ). Many cDs experiencing star formation are also experiencing 
powerful AGN outburst s. An examp l e of st ar formation in the Abell 1795 cD is 



shown in flgure [T2]from [OT)ea et al.l (|2004 l) 



There are two schools of thought regarding the origin of this star formation, 
its accompanyiiig cold gas and nebular er nission: accre tion through mergers (e.g., 
Holtzman et al.l Il992l ;i. or cooling flows (jFabianl[l993 ). Mergers are an appeal- 



ing mechanism because they are the principal route to the growth of structure 
in the context of the cold dark matter cosmogony, and they are an established 
and effective mechanism for triggering star formation in galaxies. Moreover, they 
provided a plausible alternative to classical cooling flows, that predict excessive 
;tar formation rates in cD galaxies CE^Mi). Just how kective mergers 
are at stimulating star formation in cluster cores in general and cDs in particular 
is unclear. Gas rich donor galaxies rarely inhabit the cores of clusters, perhaps 
because their gas becomes vulnerable to stripping long before they are able to 
reach the center. Furthermore, galaxies in the cores of clusters are less likely 
to host star forrn ation and nuclear activity than galaxies in other environments 
(|Best et al.ll2007l V which augurs against merger induced star formation. When 
star formation is observed in cDs, with rare exceptions, it is centered on a cooling 
flow. Most importantly, star formation rates approach or agree with the upper 
limits on cooling rates frq r n XMM-Newton and Chandra ( McNam ara et al. 20041 : 
Hicks fc Mushotzkvl[2005l : lRaffertv eralll2006l : ISalome et al.„200a ). The velocity 



structure of the cold gas fueling star forma tion is also consi s tent with having re - 
cently condensed from a static atmosphere ( Jaffe et al. 20051 : Salome et al. 20061 ). 

The emerging consistency between star formation and cooling rates is note- 
worthy in the broader context of galaxy formation. The mechanism suppressing 
cooling flows and regulating the growth of cD and gE galaxies may also be re- 
sponsible for the exponential decline in the luminosity function of bright galax- 
ies, and perhaps th e relationship between bulge luminosity and black holes (e.g., 



Benson et al. 20031 ). The possibility that these two fundamentally important 



properties of galaxies may be explained using a single mechanism has led to the 
development of new galaxy formation models that incorporate AGN feedback. 
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Cooling flows should be relatively clean examples of galaxy and supermassive 
black hole growth through the accretion of cooling gas in large halos. 

7.2 Galaxy Formation Models 

Based on the Millennium Run dark matter simulation, Croton et al. ( 20061 ) mod- 



eled the formation and evolution of galaxies and their supermassive black holes 
in a concordance CDM cosmology. Their model combined numerical and semi- 
analytic techniques, with the goal of understanding the effects of AGN feedback 
on galaxies of different masses. In this model, the AGN is powered by continu- 
ous accretion, scaling in proportion to the supermassive black hole and hot gas 
masses divided by the Hubble time. In keeping with earlier work, they follow 
galaxies through a "quasar mode" of halo merging and gas accretion that leads 
to rapid bulge and black hole growth at z = 2 - 3. This is followed by a more qui- 
escent "radio mode" that suppresses cooling flows at late times. The suppression 
becomes increasingly effective in more massive halos and is able to stop cooling 
entirely in galaxy halos with virial temperatures T > 3 x 10^ K (i.e., groups 
and clusters) from z = 1 to the present. cD galaxies continue to grow slowly at 
late times by dissipationless "dry mergers" lacking star formation. Radio mode 
feedback develops in a static hot halo driven by black hole accretion at a strongly 
sub-Eddingto i i rate that adds negligibly to the black hole mass at late times. 
Croton et all (|2006l l consider feeding the nucleus with cold cloud and Bondi ac- 



cretion from the hot halos, but lacking a feedback prescription, the outcome of 
their model is insensitive to the mode of accretion. 

Using smooth particle hydrodynamic simulations, Siiacki &: Springel ( 20061 ) 
focused on the formation and evolution of the cD and the ICM surr o undin g 
it in the presence and absence of AGN feedback. Like Croton et al. ( 20061 ). 



ISijacki k Springell (|2006l ) distinguish between major galaxy and black hole growth 
(BHAR model) at early times and more gentle AGN a ccretion and r a dio bu bble 
feedback (Magorri a n mo de) at late times. In both the Croton et al. ( 20061 ) and 
ISijacki fc Springell tood ) models, the correlation between black hole mass and 
bulge mass is imprinted during the quasar era and changes little at late times, 
when the suppression of cooling flows becomes important. This model assumes 
bubble injection repeats every 10® year with the energy per bubble scaling with 
the halo mass to the 4/3 power in the BHAR phase, and scaling with the black 
hole accretion rate during the Magorrian phase. In the absence of bubble heating, 
massive galaxies grow to unrealistically large sizes. The introduction of hot buoy- 
ant bubbles reduces cooling flows substantially, but cooling and star formation 
does not cease entirely at redshifts below z = 1 . This behavior is qualit atively 
consistent with observations of real clusters. The Sijacki &: Spriiieell (2006) model 
successfully produces entro py pedestals in the cores of clusters ( Voit &: Donahuel 
2OO5I : Donahue et al. 20061 ). However, their temperature and density profiles are 



flatter than those in real cooling flow clusters with active AGN. 

These and other studies demonstrate that periodic AGN outbursts with realis- 
tic energies and duty cycles are able to suppress cooling flows and to recover the 
observed exponential turnover in the luminosity function of bright galaxies. They 
may also play a role in the development of the black hole bulge mass relation. 
Despite these impressive results, a great deal of the picture is missing in these 
simulations, including the "microphysics" of heating and of course the "macro- 
physics" of star formation, both of which are poorly understood. These models 
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also lack a working feedback prescription that reproduces the observed levels of 
cooling, AGN activity, and star formation in gEs and cDs at late times. They 
also fail to reproduce the observed temperature and density profiles of cooling 
atmospheres. A complete model for AGN heating should include a physically 
realistic treatment of the behavior of radio jets and their interaction with the 
ICM. 

It is worth noting an important difference between the mode of energy output 
from AGN during the quasar epoch and present day cooling flows that is not 
explicitly captured in these models. Quasars are thought to be triggered by cold 
gas funneled onto the nucleus during halo mergers in the early Uni verse, lead 



ing t o the rapid buildup o f bulges and supermassive black holes (e.g., iBlandford 



19991 : ISprineel et al.l [2OO5I I . Only a small fraction of the ~ 10^^ '^^ erg s ^ of 



power radiated by this process must couple to the gas in order to drive an out- 
flow capable of suppressing star formation and black hole growth. Roughly the 
same driving power ~ 

1043-45 

erg s ^, is required to offset a massive cooling flow. 
The enthalpy and shock energy released by cooling flow AGN, even in cDs with 
relatively weak radio sources, can rival the power output of a quasar. However, 
the bulk of the power emerging from an FR I source in a cooling flow is mechan- 
ically coupled to the gas through shocks and cavity enthalpy, and is not released 
as radiation. Apparently, as the specific accretion rate throttles down from the 
Eddington regime during the quasar era to the strongly sub-Eddington regime at 
late times, the power output switches from being radiation dominated to being 
mechanically dominated. This shift i s qualitatively consist ent with the observed 
behavior of low-mass X-ray binaries ( Churazov et al. 20051 ). 



7.3 Supermassive Black Hole Growth as a Consequence of Feed- 
back 

The demanding power requirements of the long-term suppression of cluster cool- 
i ng flows by AGN shou ld produce black holes with masses exceeding 10^ Mq 
(|Fuiita fc Eeiprichll2004l V Expressing the cooling luminosity in terms of the clas- 



sical cooling rate, Lcooi = 2.5MA;T/(;umH), multiplying by the age of the cooling 
flow, t, and equating to the energy released by the black hole, rjMhC^, gives the 
minimum black hole mass required to stop cooling as 

1 5 ^2 

Mh = - — TrMt ~ 2.25 X 10"'' Mt, 

for gas cooling from 5 keV and ij = 0.1. Cooling proceeding at a rate of M = 
100 M0 year~^ since z = 1, which corresponds to a lookback time oft = 7.7 x 10^ 
year in a concordance cosmology, would form a 1.7 x 10^ M© black hole. The 
power output required to quench a 1000 Mq year~^ cooling flow (e.g., Abell 1835) 
would produce a 1.7 x 10^ black hole, rivaling the largest known black hole 
masses. Evidently, quenching a large cooling flow over the lifetime of a galaxy 
cluster could produce outsized supermassive black holes in cDs, even if jet heating 
is efficient. NGC 1275 is one of the only cluster cooling flow cDs with a measured 
black hole mass. Using molecular H2 emission line speeds, Wilman et al.l ( 20051 ) 



measured a mass of 3 x 10^ Mq. Based on the formula above, AGN feedback 
can quench a 150 - 200 Mq year~^ cooling flow without producing an excessively 
massive black hole. This may leave little room for an early buildup of the black 
hole and subsequent long-term accretion from the cooling flow. It also requires a 
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relatively high efficiency for converting accreted mass to jet power if the AGN is 
suppressing cooling. This level of growth could cause significant departures from 
the established black hole mass versus bulge mass (veloc ity dispersion) relation 
extrapolated to cD galaxy luminosities ( Lauer et al.l 2007 ). Future measurements 
of black hole ma sses in cD galaxies using very large aperture telescopes (e.g., 
Lauer et al.ll2007l ) will be a sensitive probe of the accretion and feedback history 
of clusters. 



7.4 Accretion Mechanism 

The accretion mechanism is a critical piece of any operational feedback loop. 
There are three broad categories: accretion from the hot atmosp here s urrounding 
the black hole through the Bondi mechanism (e.g., i Allen et al. 20061 ). accretion 
of col d clouds as t he result of stripping or coolin g blobs of gas from a cooling flow 
(e.g., Soker 20061 ) . and accretion of stars (e.g.. IWang &: Hu 20051 ). The actual 
mechanism must operate in the sub-Eddington accretion regime, it must release 
most of its energy by mechanical winds or outflows and not by radiation, and it 
must be responsive enough to prevent cooling and heating catastrophes. 

Bondi accretion, with a rate th at scales as neT'"^/^, is simple, natural, and 
in principle easy to regulate (e.g., IChurazov et al.ll2002l ). The local atmosphere 
responds to AGN heating by expanding, which lowers the local gas density and 
the accretion rate along with it. As the gas radiates away its energy, the atmo- 
sphere contracts and compresses, and the accretion rate rises. For 7 = 5/3, the 
Bondi accretion rate depends on gas properties only through the entropy, which 
is affected directly by heating and cooling. The accretion rate depends on the 
square of the black hole mass and properties of the atmosphere within the radius 
of influence of the black hole, all of which are difficult or impossible to measure 
with current instrumentation. Bondi accretion appears to be energetically plau- 
sible in lower luminosity gE galaxies (| Allen et al.ll2006ll. but is probably unable 
to power the largest outbursts in clusters of galaxi^TlSafel^iii). 

Ample supplies of cold gas are available in cDs to fuel the AGN. If this gas 
were supplied by stripping or mergers, it would arrive sporadically in a manner 
unrelated to the properties of the cooling atmosphere. This makes the accretion 
rate difficult to tune, unless the gas is stored in a disk, which is somehow regulated 
in the vicinity of the black hole. Gas cooling out of the IC M would na turally be 
subject to thermostatic control by periodic outbursts (e.g., Soker 20061 ). Even so, 
understanding how energy emerging from a region smaller than the size of the 
solar system is able to regulate flows on the vast scales of clusters is a monumental 
challenge. 



7.5 Observational Constraints on Feedback Models 

We pointed out earlier that the conditions in the cores of clusters lend themselves 
to detailed study of the rudiments of galaxy formation. Because the entire cycle 
of heating, cooling, feedback, and star formation can be explored there in great 
detail , cooling flows provide a standard against which models of galaxy formation 
(e.g., ISiiacki fc Soringell I2OO6I ) can be tested. These models should satisfy the 
following conditions: 

1. A small fraction, < 10%, of cDs at 2 ~ are experiencing star formation 
perhaps as a consequence of the inability of AGN to balance radiative cool- 



Heating Hot Atmospheres with AGN 



41 



3. 



ing at all times from their dense, ~ 10~^ cm~^, atmospheres with cooling 
times <C 10^ year. Infall of cold gas by mergers or stripping is another po- 
tential source of cold gas. However, it might be difficult to understand why 
stripping occurs preferentially in the cores of clusters with short cooling 
times. 

In both gEs and cDs, AGN heating scales in proportion to the cooling 
luminosity, as expected in an operational feedback loop. It is not clear 
whether or how the level of feedback scales with black hole mass or halo 
mass. 



Bondi accretion may power the AGN in gEs ( Allen et al.l 20061 ) but is prob- 
ably unable to do so in more massive cD galaxies. Only the most powerful 
cluster outbursts require accretion rates approaching ~ 1 Mq year~^; typi- 
cal rates are ~ 10-2 Mq year , which is a small fraction of the Eddington 
accretion rate of a 10^ Mrr^ black hole. 



4. 



The je t model most directly associated with the M agorrian mode of lSiiacki Sz Springe! 
(|2006l l and the radio mode of lCroton et al.l (|2006l l must account for the enor- 
mous range in radiative efficiency and large jet kinetic energy observed in 
cDs. It should also recover the observed temperature, density, and en tropy 
profiles in the cores of clusters (iBorgani I2OO4I : IVoit fc Donahu3l200,'j ). 

The existence of very powerful AGN outbursts and the persistent energy 
demands of cooling imply substantial black hole growth at late times in 
cD galaxies. There is tantalizing new evidence that cDs and their central 
black holes grow at an average rate that cr udely follows the sl ope of the 
black hole mass versus bulge m ass relation ( Haring &: Rix 20041 ) once star 
formation is taken into account ( Raffertv et al. 20061 ). However, bulges and 
black holes do not always grow in lock step. Other lines of evidence suggest 
that the slope flattens in cD galaxies ( Lauer et ahlboo"?! ). 

Star formation parameterized with a Schmidt- Kennicutt law is probably a 



reaso nable approximation in these systems (e.g.. lEgami et al.ll2006l : lMcNamara et al. 



20061 ). although this issue is in need of further study. Disk formation, which 



is a s taple of semianalytic galaxy formation rn odels, is rare or short-lived 
(e.g., iHeckman etla]ll989l : ISalome et al.l l2006l ;i. Energy pumped into the 
hot gas by st arburst winds, 5, 10^^ erg s~ \ is negligible on the scale of the 
cooling flow (jMcNamara et al.ll20oH2006l ). but may be important near the 
nucleus of the cD where fueling of the AGN is actively occurring. 

7. Finally, the observed level and spatial distribution of chemical enrichment, 
and abundance ratios in the hot gas must constrain the history of star 
formation in the cores of clusters and the level of mixing generated by 
AGN outbursts and mergers. 



8 CONCLUDING REMARKS 



The editors of Discover Magazine's Winter 2007 special issue "Unseen Universe" 
posed the following question to Martin Rees: "Is there any particular image that 
you saw recently that reminded you anew of just how much progress we've made?" 
Rees's response provides a succinct history of this field. "I've been especially im- 
pressed by the X-ray images of galaxy clusters which are now becoming available 
from the Chandra satellite and other instruments. We see gas being churned up 
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by explosions and huge black holes in the center of the cluster. We see how it's 
cooling down and how the cooling is being balanced by tremendous outbursts of 
jets and bubbles of hot gas. This is something that most people didn't suspect 
was happening until these images revealed it." 

If many people did not suspect this was happening before Chandra revealed 
it, why not? The Rosat observatory had alread y established that radio sources are 
inter acting with the hot gas surrounding them (jBohringer et al lll99.4l ICarilli et al 



199411. and it had been pointed out previously t hat AGN (e .g., 'Pedlar et al."l99C 



Baum fc 0'Dealll99ll : [Owen fc E ilek 1998; Tuc ker &: Davidl[i9 97: Binney & Tabo] 
19951 1 and thermal conduction (jRosner &: Tuckeiill989l : Meiksinill988l ) could offset 



cooling in some clusters. In our view, these solutions were not widely embraced 
because, like the classical cooling flow model itself, they lacked persuasive obser- 
vational and theoretical support. As is often true in astronomy, the impasse was 
broken by new, high-resolution instrumentation. XMM-Newton spectroscopy re- 
vealed that hot gas in clusters cannot be cooling at the classical rates. At about 
the same time, high resolution Chandra images identified AGN feedback as the 
probable heating agent. But the real situation is more complicated. The data 
show that AGN heating is more subtle than early theoretical studies envisioned. 
Heating is not a violent, local process. It is gentle and spatially dispersed. There 
is scant evidence for constant density cores or central temperature spikes, and 
entropy inversions as some nuclear heating models predicted. To our knowledge, 
no one anticipated cool rims surrounding the cavities, ripples, ghost cavities, and 
quasar-like outbursts that are barely audible in the radio and optically faint in 
the nucleus. The data were inadequate and so was the physics. 

Much remains to be done. We do not understand how jets are powered and 
what they are made of, how the putative feedback loop works, how efficiently jets 
heat the gas, and we still cannot be sure that they are the sole heating agent or 
even the principal one. These issues will take time to resolve, but a great deal 
of progress has already been made, notably in ever more realistic jet simulations 
(e.g.. Figure [TT]) that can now be tested against high quality data. They must 
be resolved in order to understand the role of supermassive black holes in galaxy 
and large-scale structure formation. We conclude by listing several issues and 
avenues of research that we believe will lead to substantial progress in this field. 



FUTURE ISSUES 



1. XMM-Newton spectroscopy has excluded mass deposition at the classical 
rates in cooling flows. However, deposition at substantial levels comparable 
to observed star-formation rates have not been ruled out. A search for 
Fe XVII and other lines at levels that are consistent with observed star 
formation rates is within the grasp of XMM-Newton's reflection grating 
spectrometer, and will be within easy reach of a future Constellation-X. 
Combined with increasingly accurate star formation rates (item 5), this 
would constitute a strong test of new, feedback-based cooling flow models. 

2. It is difficult or impossible to detect cavities and shock fronts in distant 
clusters and groups beyond z ~ 0.5. In principle, radio observations, espe- 
cially at low frequencies, can probe the history of feedback and heating at 
earlier times. However, the large variation in radio synchrotron efficiency 
must first be understood and calibrated. 
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3. The time distribution of AGN jet power in galaxies and clusters, which is 
an essential part of a feedback model, is poorly understood. Existing X-ray 
cavity surveys are limited in size and suffer strong selection biases. A large, 
unbiased search for cavities and shocks in a flux or volume limited sample 
of groups and clusters is required to more accurately determine the average 
AGN heating rate. 

4. The environmental factors that trigger AGN outbursts, particularly the 
most powerful ones in cD galaxies (> 10^^ erg), are poorly understood. 
High resolution imaging and spectroscopy of the stars and gas in nuclei of 
galaxies, made at a variety of wavelengths, will characterize the properties 
of nuclear gas disks, stellar cusps, and cores. This knowledge should provide 
a deeper understanding of the accretion and feedback process, and hopefully 
new insight as to why many systems are as powerful as quasars and yet they 
don't reveal themselves as such. Most importantly, measurements of black 
hole masses in cD galaxies using future large aperture telescopes will place 
restrictive limits on the history of AGN feedback in clusters. 

5. Limits on the rate of cooling in the cores of galaxies and clusters have be- 
come tight enough to warrant better measurements of star formation rates 
and histories. This will require precision photometry over a broad pass- 
band (UV to far IR) and careful acc ounting for dust and positive feedb ack 
from radio triggered star formation ( De Youn3 19951 : Q'Dea et al. 2004 : see 



Figure [T2]) . Combining this information with black hole growth rates esti- 
mated from cavities and shock fronts will yield new insight into how bulges 
and supermassive black holes grow at late times. 



6. New high fidelity jet models (e.g., iHeinz et al.l l2006l see Figure [TT]) com- 



bined with sensitive X-ray measurements of cavities and shock fronts will 
increase our understanding of jet dynamics, energetics, content, radiative 
efficiencies, and ultimately how jets form. The form of energy deposition 
by jets is a critical component of any feedback model. 

7. Models are needed for the fueling and triggering of AGN outbursts, includ- 
ing the part played by mergers. 

8. How cavity enthalpy and weak shock energy is dissipated, how efficiently it 
heats the gas, and where the heat is deposited are fundamentally important 
questions that have not been satisfactorily answered. 

9. The contribution of AGN outbursts to deviations from the expected scaling 
relati onships between mass, temperature, and X-ray luminosity of clusters 
fe.g.. lO'Hara et al.l[2006l : iGitti et aDl2007l ^ is poorly understood. The en- 



ergy required to quench a large cooling flow over its lifetime is comparable to 
the 1 keV per baryon necessary to preheat a cluster. AGN may contribute 
significantly to preheating. 
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Figure 1: Left. Chandra X-ray ima ge of M87 showing structure in the hot gas 
associated with the AGN outburst ( Forman et al. 2003). Several smah cavities 



and filaments of uplifted gas are visible. X-ray emission from the radio jet is 
seen at the center. Right. Chandra X-ray image of the Perseus cluster exposed 
for 900,000 s from Fabian et al. (2006). Two inner cavities containing the active 
radio lobes and two outer ghost cavities are seen. The sound waves (weak shocks) 
are visible as a series of circular ripples surrounding the cavities. The central blue 
structure is absorption by a foreground galaxy. 
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Figure 2: Hubble Space Telescope visual image of the MS0735. 6+7421 cluster 
superposed with the Chandra X-ray image in blue and a radio image from the 
Very Large Array at a frequency of 330 MHz in red. The X-ray image shows an 
enormous pair of cavities, each roughly 200 kpc in diameter that are filled with 
radio emission. The radio jets have been inflating the cavities for 10^ years with 
an average power of < 2x 10^^ erg s~^. The displaced gas mass is < 10^^ Mq. The 
cavities and radio source are bounded by a weak shock front. The cavities are well 
outside the central galaxy and cooling region of the cluster. The supermassive 
black hole grew by at least < 3 x 10® M0 during the outburst. 
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Figure 3: Distribution of cavity distance from the nucleus of the host cD galaxy. 
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Figure 4: Distribution of cavity sizes. An ellipse was fitted to each cavity and its 
is size expressed as Vab, where a and b are the semiaxes. 
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Figure 5: Distribution of the ratio of cavity size to nuclear distance. See Figures 
3 and 4. 
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Figure 6: Distribution of cavity buoyancy ages. 



Heating Hot Atmospheres with AGN 57 



10^ F 




Figure 7: Total radio luminosity (10 MHz - 10 GHz) plotted against jet power 
(4py/tbouy) taken from Laura Bizan's PhD thesis (2007). The diagonal Unes 
represent ratios of constant jet power to radio synchrotron power. Jet power 
correlates with synchrotron power but with a large scatter in their ratio. Radio 
sources in cooling flows are dominated by mechanical power. The radio measure- 
ments were made with the Very Large Array telescope. 
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Figure 8: Cavity power of the central AGN plotted against the X-ray luminosity of 
the intraclu ster medium (ICM) w ithin the cooling radius, after correcting for mass 
deposition ( RafFertv et al. 20061 ). The symbols and wide error bars denote values 
of cavity power calculated using the buoyancy timescale. Short and medium 
width error bars denote the limits of the cavity power calculated using the sound 
speed and refill timescales, respectively. Diagonal lines denote equality between 
heating and cooling rates assuming pV, 'ipV, and 16pV of energy per cavity, 
respectively. 
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Figure 9: Cavity po wer versus cooling power for nearby giant elliptical galaxies 
(jNulsen et alJ 120071 ^. Cooling power is the X-ray luminosity from within the 
projected radius where the cooling time is 7.7 x 10^ years. Cavity powers are 
determined using an energy of pV per cavity and a range of cavity age estimates 
(see Figure 8). The dashed lines show equality for cavity energies of "pV ^ ApV, 
and 16pV, top to bottom. All but one system lie above the 4:pV line, indicating 
that radiative losses can bebalanced by AGN power. 




Figure 10: Buoyant cavity. As the cavity rises, gas falls inward to fill the space, 
turning gravitational potential energy into kinetic energy in the cavity's wake. 
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Figure 11: S i mulat ion of jets interacting with the intracluster medium (ICM) from 



Heinz et al 



( 20061 ) . Radio synchrotron emission is in blue and X-ray emission in 
red. This simulation produces a realistic cocoon shock, cavities, and ripples. The 
strong gas feature to the lower right is an unrelated merger shock. Motion of 
ICM owing to ongoing cluster growth and wobbling of the Cygnus A-like jet both 
have a significant impact on the outcome of the model. 
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Figure 12: Hubble Space Telescop e images of the ceii tral 19 x 19 arcsec (24 x 24 
kpc) of the cD in Abell 1795 from lO 'Pea et al.l (|2004l ): Left: Far-UV continuum, 
middle: Ly emission with radio contours superposed, and right: I-band contin- 
uum. Note the bright knots of star formation and Ly emission that are seen 
preferentially along the radio lobes. A dust lane is evident in the I-band image. 
The star formation rate in this system is 10 - 30 M© year~^. 



